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1  SUMMARY 


Recent  research  has  shown  that  one  possible  explanation  to  bypass  transition  is  the  transient  growth 
theory.  However,  there  has  not  been  any  direct  numerical  simulation  study  on  transient  growth  in 
hypersonic  boundary  layers.  Furthermore,  it  is  not  known  how  the  optimal  disturbances  computed 
by  the  transient  growth  theory  are  generated  by  surface  roughness.  It  is  also  not  clear  what  the  role 
of  freestream  disturbances  are  in  the  transient  growth  theory.  These  are  important  issues  related  to 
the  receptivity  of  transient  growth,  which  need  to  be  resolved.  The  understanding  of  such  a  complex 
flow  phenomena  can  best  be  obtained  through  a  coordinated  study  involving  both  computational 
and  theoretical  analyses.  Direct  numerical  simulation  (DNS)  is  uniquely  suited  for  such  studies 
because  the  complete  receptivity  process  and  transient  growth  can  be  simulated  with  minimum 
simplification  assumptions.  Supported  by  AFOSR  under  the  Grant  number  FA9550-04- 1-0029  from 
January  2004  to  December  2006,  our  group  at  UCLA  has  demonstrated  that  our  unique  high-order 
shock  fitting  simulation  approach  is  a  powerful  tool  in  studying  supersonic  and  hypersonic 
boundary-layer  stability  and  transition  physics.  We  have  conducted  extensive  DNS  studies  on  the 
flow  mechanisms  of  hypersonic  boundary  layer  receptivity  and  stability.  Our  computer  code  is  fully 
validated  by  comparing  numerical  simulation  solutions  with  available  experimental  or  theoretical 
results  on  hypersonic  boundary  layer  receptivity  and  stability. 

The  objective  of  the  three-year  research  is  to  conduct  DNS  studies  of  hypersonic  boundary  layer 
receptivity,  transient  growth  and  transition  with  surface  roughness.  The  main  approach  of  the 
proposed  research  is  to  use  DNS  as  a  research  tool  to  study  the  boundary  layer  receptivity  and 
transient-growth  mechanisms  in  hypersonic  flows,  including  the  development  of  numerical 
algorithms  and  parallel  computer  codes  of  higher  order  numerical  methods  for  the  simulation  of  hypersonic 
flows  with  surface  roughness  of  finite  heights.  During  the  three-year  period,  we  have  conducted  DNS 
studies  on  the  hypersonic  boundary  layer  flows  over  fiat  plates  and  blunt  cones.  A  new  high-order 
cut-cell  method  has  been  developed  for  the  numerical  simulation  of  hypersonic  boundary  layer 
transition  with  finite  height  surface  roughness.  The  method  has  been  applied  to  the  numerical 
simulations  of  two-dimensional  hypersonic  flows  over  a  flat  plate.  Furthermore,  the  stabilization 
effect  of  the  surface  porous  coating  over  a  fiat  plate  is  extensively  studied  by  series  of  numerical 
simulations.  We  also  collaborate  with  Prof  Tumin  in  the  University  of  Arizona  to  compare 
numerical  and  theoretical  results  on  receptivity  of  a  Mach  5.92  flow  over  a  fiat  plate  to  wall 
blowing-suction,  and  to  deeply  analyze  the  nonparallel  flow  effects.  Our  numerical  studies  have 
been  validated  to  be  of  high  accuracy  and  led  to  further  understanding  of  hypersonic  boundary  layer 
receptivity  mechanism.  Such  understanding  can  lead  to  better  tool  for  the  prediction  and  control  of 
high  speed  boundary  layer  transition. 

The  main  research  contributions  are: 

1.  We  have  developed  a  new  high-order  cut-cell  method  for  the  numerical  simulation  of 
hypersonic  boundary  layer  transition  with  finite  height  surface  roughness.  Parallel  computer  codes 
of  the  numerical  method  have  been  implemented.  The  new  high-order  cut-cell  method  has  been 
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tested  in  the  computations  of  several  one  and  two-dimensional  hyperbolic  and  elliptic  equation 
in  irregular  domains.  The  results  show  that  up  to  fourth-order  accuracy  in  bothZ^  and  norm 
can  be  obtained  for  the  current  cut-cell  method  for  both  problems. 

2.  We  subsequently  have  applied  a  third-order  cut-cell  method  to  the  two-dimensional 
compressible  Navier-Stokes  equations  for  simulating  roughness  induced  receptivity  for 
hypersonic  flow  over  a  flat  plate  with  a  blowing  and  suction  slot  placed  near  the  leading  edge. 
The  surface  roughness  height  is  approximately  half  of  the  local  boundary  thickness.  The  results 
show  that  the  roughness  location  plays  an  important  role  in  the  developments  of  mode  S  initially 
excited  by  the  blowing-suction  actuator.  Mode  S  is  destabilized  only  when  the  roughness 
element  is  located  close  upstream  of  the  synchronization  point.  On  the  other  hand,  when  the 
roughness  element  is  downstream  of  the  synchronization  point,  mode  S  is  stabilized.  This 
happens  even  when  the  roughness  is  still  within  the  unstable  region  of  mode  S.  Therefore,  the 
synchronization  point  and  finite  roughness  location  is  critical  to  the  receptivity  process. 

3.  We  have  extensively  studied  the  stabilization  of  a  Mach  5.92  fiat-plate  boundary  layer  using 
local  sections  of  porous  coating  by  series  of  numerical  simulations.  The  stabilization  effect  of 
porous  coating  on  hypersonic  boundary  layers  over  flat  plates  and  cones  has  been  successfully 
demonstrated  by  theoretical  analyses,  experiments,  and  numerical  simulations.  It  has  been  found 
that  porous  coating  slightly  destabilizes  the  Mack  first  mode  whereas  it  significantly  stabilizes 
the  Mack  second  mode.  In  previous  studies,  porous  coating  covers  either  the  entire  fiat  plate  or 
the  surface  around  half  the  cone  circumference.  The  effect  of  porous  coating  locations  on 
boundary  layer  stabilizations  has  not  been  considered.  Furthermore,  the  destabilization  effect  of 
porous  coating  on  the  Mack  first  mode  has  not  been  studied  in  detail.  In  our  numerical 
simulations,  we  focus  on  the  effect  of  porous  coating  locations  and  the  destabilization  of  the 
Mack  first  mode.  It  is  found  that  the  synchronization  point  of  mode  S  and  mode  F  plays  an 
important  role  on  boundary  layer  stabilization.  Disturbances  are  destabilized  when  porous 
coating  is  located  upstream  of  the  synchronization  point  whereas  they  are  stabilized  when 
porous  coating  is  downstream  of  the  synchronization  point.  For  felt-metal  porous  coating, 
destabilization  of  the  Mack  first  mode  is  significant. 

4.  We  collaborate  with  Prof  Tumin  in  the  University  of  Arizona  to  compare  numerical  and 
theoretical  results  on  receptivity  of  the  Mach  5.92  flow  over  a  flat  plate  to  wall  blowing-suction, 
and  to  deeply  analyze  the  nonparallel  flow  effects.  The  perturbation  flow  field  is  decomposed 
into  normal  modes  with  the  help  of  the  multimode  decomposition  technique  based  on  the  spatial 
biorthogonal  eigenfunction  system.  The  decomposition  allows  filtering  out  the  stable  and 
unstable  modes  hidden  behind  perturbations  having  another  physical  nature.  The  results  indicate 
that  nonparallel  flow  effects  stabilize  the  second  mode  whereas  destabilize  the  first  mode. 

5.  We  have  also  conducted  numerical  simulation  studies  of  Stetson’s  Mach  5.5  experiments  on  the 
nose  bluntness  effects.  Three  nose  radii  of  0.0156in,  0.5in,  and  1.5in  are  used  to  study  the  nose 
bluntness  effects  on  boundary  layer  instability  and  transition.  The  simulation  results  show  that 
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the  increase  of  nose  bluntness  leads  to  the  substantial  decrease  in  local  Reynolds  numbers  along 
the  edge  of  the  boundary  layers. 

2  RESEARCH  OBJECTIVES 


The  prediction  of  laminar-turbulent  transition  in  supersonic  and  hypersonic  boundary  layers  is  a 
critical  part  of  the  aerodynamics  and  heating  analyses  for  the  development  of  hypersonic 
transportation  vehicles  and  re-entry  vehicles.  The  key  for  an  accurate  transition  prediction  is  the 
understanding  of  the  physical  mechanisms  that  lead  to  transition.  However,  many  important 
physical  mechanisms  leading  to  hypersonic  boundary  layer  transition  are  currently  still  not  well 
understood.  Recent  research  has  shown  that  one  possible  explanation  to  bypass  transition  is  the 
transient  growth  theory.  However,  there  has  not  been  any  direct  numerical  simulation  study  on 
transient  growth  in  hypersonic  boundary  layers.  Furthermore,  it  is  not  known  how  the  optimal 
disturbances  computed  by  the  transient  growth  theory  are  generated  by  surface  roughness.  It  is  also 
not  clear  what  the  role  of  freestream  disturbances  are  in  the  transient  growth  theory.  These  are 
important  issues  related  to  the  receptivity  of  transient  growth,  which  need  to  be  resolved.  Due  to  the 
difficulties  in  conducting  hypervelocity  experiments  and  the  complexity  of  hypersonic  flows, 
fundamental  hypersonic  studies  increasingly  rely  on  the  use  of  DNS  as  a  research  tool.  In  recent 
years,  DNS  has  become  a  powerful  tool  in  the  studies  of  the  stability  and  receptivity  of  supersonic 
and  hypersonic  boundary  layers. 

Supported  by  AFOSR,  we  have  developed  high-order  shock-fitting  DNS  methods  [1],  which  can  be 
directly  applied  to  hypersonic  boundary  layers  over  realistic  blunt  bodies  with  the  effects  of  nose 
bluntness,  the  presence  of  bow  shock  waves,  and  the  real-gas  effects  at  high  temperatures.  In  the 
past  several  years,  we  have  studied  the  receptivity  and  stability  of  a  number  of  2-D  and  3-D 
hypersonic  flows  over  flat  plates  [2-6],  and  blunt  cones  [7-10]. 

The  main  objectives  of  this  research  are  to  conduct  extensive  direct  numerical  simulation  studies  on 
the  early  stage  of  hypersonic  boundary  layer  receptivity,  transient  growth,  and  transition  for  blunt 
cones  and  flat  plates,  and  to  develop  numerical  algorithms  and  parallel  computer  codes  of  higher 
order  numerical  methods  for  the  simulation  of  hypersonic  flows  with  surface  roughness  of  finite 
heights.  The  effects  of  surface  roughness,  nose  bluntness,  and  freestream  disturbances  are  studied 
by  numerical  simulations. 

The  project  has  been  carried  out  in  three  years  with  the  following  research  tasks: 

1.  Development  of  new  high-order  direct  numerical  simulation  methods  and  3-D 
computer  codes  to  simulate  arbitrary  surface  roughness  effects.  A  new  high-order 
cut-cell  method  for  the  numerical  simulation  of  hypersonic  boundary  layer  transition 
with  finite  height  surface  roughness  [3,  4].  Parallel  computer  codes  of  the  numerical 
method  have  been  implemented.  The  new  high-order  cut-cell  method  has  been  tested  in 
the  computations  of  several  one  and  two-dimensional  hyperbolic  and  elliptic  equation 
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in  irregular  domains.  The  results  show  that  up  to  fourth-order  aeeuraey  in  bothT^  and 
norm  ean  be  obtained  for  the  eurrent  eut-eell  method  for  both  problems.  We 
subsequently  have  applied  a  third-order  eut-eell  method  to  the  two-dimensional 
compressible  Navier-Stokes  equations  for  simulating  roughness  induced  receptivity  for 
hypersonic  flow  over  a  flat  plate  with  a  blowing  and  suction  slot  placed  near  the 
leading  edge.  The  surface  roughness  height  is  approximately  half  of  the  local  boundary 
thickness. 

2.  DNS  studies  of  the  stabilization  of  a  Mach  5.92  flat-plate  boundary  layer  using 
local  sections  of  porous  coating.  In  our  numerical  simulations  [2],  we  focus  on  the 
effect  of  porous  coating  locations  and  the  destabilization  of  the  Mack  first  mode.  It  is 
found  that  the  synchronization  point  of  mode  S  and  mode  F  plays  an  important  role  on 
boundary  layer  stabilization.  Disturbances  are  destabilized  when  porous  coating  is 
located  upstream  of  the  synchronization  point  whereas  they  are  stabilized  when  porous 
coating  is  downstream  of  the  synchronization  point.  For  felt-metal  porous  coating, 
destabilization  of  the  Mack  first  mode  is  significant. 

3.  Collaborate  with  Prof.  Tumin  in  the  University  of  Arizona  to  compare  numerical 
and  theoretical  results  on  receptivity  of  the  Mach  5.92  flow  over  a  flat  plate  to  wall 
blowing-suction,  and  to  deeply  analyze  the  nonparallel  flow  effects  [6].  The 

perturbation  flow  field  is  decomposed  into  normal  modes  with  the  help  of  the 
multimode  decomposition  technique  based  on  the  spatial  biorthogonal  eigenfunction 
system.  The  decomposition  allows  fdtering  out  the  stable  and  unstable  modes  hidden 
behind  perturbations  having  another  physical  nature.  The  results  indicate  that 
nonparallel  flow  effects  stabilize  the  second  mode  whereas  destabilize  the  first  mode. 

4.  DNS  studies  of  nose-bluntness  and  entropy-layer  effects  on  hypersonic  boundary 
layer  in  the  context  of  transient  growth  with  or  without  surface  roughness.  We 

conduct  DNS  investigations  of  transition  reversal  due  to  the  nose  bluntness  effects  for 
hypersonic  boundary  layer  over  blunt  cones  and  blunt  fiat  plates  [10,  11].  A  large 
number  of  test  cases  of  various  nose  bluntness,  Mach  numbers,  surface  roughness,  and 
Reynolds  numbers  are  studied.  We  study  the  nose  bluntness  effects  in  the  context  of 
receptivity,  transient  growth,  and  subsequent  development  of  instability  waves. 


3  RESEARCH  ACCOMPLISHMENTS 


The  research  project  described  in  the  previous  section  was  carried  out  in  a  three  period  from 
January  1,  2007  to  December  31,  2009.  The  three-year  research  project  supported  by  this  grant  has 
led  to  the  publications  of 

1 .  Six  papers  in  the  archive  journals  or  book  chapters. 
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2.  Eighteen  conference  papers  (mainly  AIAA  papers), 


More  importantly,  the  research  has  produced  a  large  amount  of  new  results,  and  has  led  to  new 
understanding  of  a  number  of  hypersonic  boundary  layer  receptivity  mechanisms.  The 
understanding  of  the  receptivity  process  is  the  key  to  improve  current  hypersonic  boundary  layer 
transition  prediction  methods.  In  the  following  sections,  the  major  research  accomplishments  are 
summarized,  followed  by  a  more  details  discussions  of  some  of  our  new  results  in  hypersonic 
boundary  layer  receptivity. 

3.1  Development  of  a  new  high-order  cut-cell  method 
Publications:  [1],  [7],  [18] 

During  the  three-year  project,  we  developed  a  new  high-order  cut  cell  method  for  computing 
boundary  layer  transition  with  surface  roughness.  Due  to  the  difficulty  in  grid  generation  for 
numerical  simulation  of  high-speed  flow  with  arbitrary  surface  roughness  elements,  it  is 
advantageous  to  use  a  fixed-grid  cut-cell  method  to  compute  such  flow.  The  grid  can  be  a  smooth 
curvilinear  body-fitted  one  along  a  baseline  smooth  surface  without  the  roughness.  The  actual 
surface  with  roughness  will  cut  across  the  grid  lines,  which  create  irregular  grid  points.  Most  of  the 
previous  Cartesian-grid  cut-cell  (or  sharp-interface)  methods  with  different  boundary  treatments, 
however,  are  only  first  or  second  order  accurate  at  the  interface.  The  relative  low  accuracy  at  the 
interface  may  not  be  sufficient  for  numerically  simulating  laminar-turbulent  boundary-layer 
transition  involving  surface  roughness.  For  a  simulation  study  of  hypersonic  boundary  layer 
transition  shown  in  Fig.  1,  it  is  necessary  to  consider  the  effects  of  both  a  bow  shock  and  an 
arbitrary  roughness  element  on  the  lower  surface.  More  details  can  be  found  in  [1].  Numerical 
methods  for  such  studies  need  to  be  high-order  accurate  both  in  resolving  the  complex 
shock/disturbance  interaction  and  in  resolving  a  wide  range  of  flow  time  and  length  scales.  In 
addition,  high-order  accuracy  is  required  for  the  local  finite  difference  schemes  for  local  region 
adjacent  to  the  surface  of  the  roughness  element.  In  this  paper,  we  present  a  new  high-order  cut-cell 
method  for  numerical  simulation  of  roughness  induced  transition  in  viscous  hypersonic 
boundary-layer  flows.  This  method  is  mainly  based  on  finite  difference  method,  and  is  different  in 
implementation  from  all  the  methods  discussed  above.  To  overcome  the  small  cell  problem,  we 
adopt  Gibou  and  Fedkiw’s  shifting-interpolation-stencil  strategy,  i.e.,  the  “drop  point”  approach. 
The  new  cut-cell  method  is  uniformly  high-order  for  both  regular  and  irregular  grid  points.  The 
computational  accuracy  of  our  third  and  fourth  order  accurate  cut-cell  methods  is  tested  for  the 
computations  of  Poisson  and  hyperbolic  equation  with  Dirichlet  boundary  conditions  in  irregular 
domains.  We  then  test  a  third-order  accurate  cut-cell  method  for  the  two-dimensional  Navier-Stokes 
equations.  The  new  scheme  is  used  to  simulate  steady  and  unsteady  hypersonic  boundary  layer 
flows  over  a  fiat  plate  with  an  isolated  surface  roughness  element  (Fig.  1).  The  bow  shock  generated 
from  the  leading  edge  of  the  flat  plate  is  treated  as  a  computational  boundary  and  discretized  based 
on  Zhong’s  [1]  high-order  finite-difference  shock- fitting  method.  The  unsteady  flow  over  fiat  plate 
is  excited  by  a  narrow  blowing  and  suction  slot  composed  of  15  different  frequencies  mounted  near 
the  leading  edge.  Several  multi-dimensional  cases  are  tested  to  demonstrate  the  order  of  accuracy 
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for  these  numerical  methods. 


Ma>1 


Fig.  1 .  An  example  hypersonic  flow  over  flat  plate  with  surface  roughness  induced  boundary  layer  transition. 


Figure  2  shows  the  numerical  solution  of  an  elliptic  equation  in  an  irregular  domain  given  by 
u^+Uyy  =  J^2<J(x-X(5'))J(y-7(5))(is' .  The  interfaceF  is  a  circle  defined  byx^+y^  =1/4.  The 

computational  domain  is  (x,y)  e  [-1,1]^  .  A  fourth  order  cut  cell  method  is  used  to  compute  the 

two-dimensional  example  and  the  exact  solution  for  the  case  of  80x80  grids.  This  figure  shows  a 
very  good  agreement  between  the  exact  and  numerical  solutions.  There  is  no  spurious  oscillation  for 
the  numerical  solutions  in  the  region  adjacent  to  the  interface. 


Fig.  2.  Comparison  of  the  solution,  u(x,y),  computed  by  Version  2  of  the  fourth-order  cut-cell  method 

)  for  the  two-dimensional  example  and  the  exact  )  (grid:  80x80 ):  1)  and  2)  are  contours  of  the 
solution,  3)  is  distribution  along  the  x  direction. 
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Table  1  shows  the  eomputational  errors,  in  and  D  norms  respeetively,  with  four  different  grid 
sets.  As  expeeted,  table  1  shows  that  the  orders  in  infinity  norm  are  lower  than  those  of  the  seeond 
norm.  Overall,  the  fourth-order  eut-eell  method  with  third-order  boundary  elosure  sehemes  produees 
a  fourth-order  global  aceuracy  in  U  norm  for  the  solutions  to  the  two-dimensional  problem. 

Table  1.  Computational  errors  \nU  norm  of  the  two  versions  of  the  cut-cell  methods  (O(h^)  and  OQi'^)) 
for  the  two-dimensional  test  case  of  LeVeque  and  Li  [12]. 


Cut  Cell  Method  1 

Cut  Cell  Method  2 

N 

0 

II F  II 
|pv|l2 

Ratio 

P 

0 

II F  II 

IPvllj 

Ratio 

P 

20 

0.5 

4.791  (-4) 

1.0 

1.967(-4) 

40 

0.5 

4.040(-5) 

11.9 

3.56 

1.0 

1.437(-5) 

13.7 

3.77 

80 

0.5 

2.919(-6) 

9.54 

3.25 

1.0 

9.049(-7) 

15.9 

3.99 

160 

0.5 

2.945(-7) 

9.91 

3.30 

1.0 

5.773(-8) 

15.7 

3.97 

3.2  Studied  the  stability  of  a  mach  5,92  flat-plate  boundary  layer  with  finite  height  roughness 
in  different  locations 

Publications:  [7],  [14] 

The  reeeptivity  of  a  Maeh  5.92  boundary-layer  flow  over  a  flat  plate  with  finite  surfaee  roughness  to 
periodie  two-dimensional  wall  perturbations  is  investigated  by  numerieal  simulations  and  linear 
stability  theory  (LST).  The  flow  eonditions  are  the  same  as  those  of  Maslov  et  al.’s  leading-edge 
experiment  [13].  The  steady  base  flow  is  simulated  by  solving  two-dimensional  eompressible 
Navier-Stokes  equations  with  a  eombination  of  a  fifth-order  shoek-fitting  method  and  a  third-order 
eut  eell  seheme.  Aeeuraey  of  the  numerieal  steady  base  flow  is  validated  by  eomparing  with  the 
theoretieal  self-similar  boundary-layer  solution  and  the  published  experimental  results.  The 
eharaeteristies  of  boundary-layer  wave  modes  are  identified  and  evaluated  by  eomparing  the  results 
of  LST  and  numerieal  simulations.  In  reeeptivity  simulations,  blowing-suetion  is  used  as  the 
periodie  two-dimensional  wall  perturbation,  whieh  is  introdueed  to  the  steady  base  flow  by  a 
foreing  slot  loeated  on  the  fiat  plate.  The  effeets  of  finite  surfaee  roughness  on  reeeptivity  process 
are  studied  by  considering  four  cases  of  different  roughness  location.  Fig.  3  compare  normalized 
Mach  number  and  dimensionless  streamwise  velocity  distributions  across  the  boundary  layer  at 
three  different  locations  of  x  =  0.00096  m,  0.00121  m,  and  0.000138  m.  The  numerical  results  agree 
well  with  the  experimental  results  and  the  boundary-layer  solution  near  the  plate.  However,  in  the 
region  of  q  >5,  the  numerical  results  have  a  better  agreement  with  the  experimental  results.  The 
difference  between  the  numerical  results  and  the  boundary-layer  solution  is  mainly  caused  by  the 


existence  of  the  bow  shock,  because  the  effect  of  the  bow  shock  is  neglected  in  the  calculation  of 
the  compressible  boundary-layer  equations. 


Fig.  3  Numerical  simulation  of  a  Mach  5.92  flat-plate  boundary  layer  to  two-dimensional  wall  perturbations  and 
comparisons  of  normalized  Mach  number  and  dimensionless  streamwise  velocity. 


Four  cases  are  simulated  by  using  the  high-order  cut  cell  method.  Specifically,  roughness  locations 
with  respect  to  the  synchronization  point  at  x  =  0.33  Im  are  as  follows, 

•  Case  1:  Upstream  of  the  synchronization  point  x  =  0.185m;  The  roughness  height 

=  =  0.00081m. 

2 

•  Case  2:  Upstream  of  the  synchronization  point  x  =  0.260m;  The  roughness  height 

=  =  0.001 11m. 

2 

•  Case  3:  Synchronization  point  x  =  0.331m;  The  roughness  height  =  0.00141m 

•  Case  4:  Downstream  of  the  synchronization  point  x  =  0.410m;  The  roughness  height 

/2  =  i^  =  0.00171m. 

2 


Fig.  4  Numerical  simulation  of  receptivity  of  the  Mach  5.92  flow  over  a  flat  plate  to  two-dimensional  wall 
perturbations  with  surface  roughness:  1)  amplitude  comparison  of  pressure,  2)  &  3)  comparison  of  numerical 

and  LST  growth  rates. 
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Fig.  4  shows  the  pressure  amplitude  on  the  flat  plate  surfaee  with  roughness  loeated  in  different 
loeations.  It  shows  that  for  all  four  eases,  the  development  of  two-dimensional  wall  perturbations 
eventually  result  in  the  same  type  of  instability  wave  (mode  S)  in  the  boundary  layer.  When  the 
surfaee  roughness  loeated  in  the  far  upstream  of  the  synehronization  point  (Case  1),  the  instability 
wave  develops  in  the  same  way  as  the  ease  without  roughness.  But  as  the  roughness  moves 
gradually  into  the  downstream  (Case  2),  the  growth  rate  of  the  mode  S  inereases,  whieh  refleets  ehat 
the  amplitude  of  pressure  in  the  flat  plate  surfaee  is  amplified.  Mode  S  is  destabilized  when  the 
roughness  element  is  plaeed  suffieient  elose  to  the  synehronization  point.  In  Case  3  and  4,  when  the 
roughness  loeated  in  downstream  of  synehronization  point,  the  flnite  surfaee  roughness  stabilizes 
the  mode  S  inside  the  boundary  layer.  The  amplitude  of  pressure  perturbation  is  redueed  in  further 
downstream. 

3,3  Studied  the  stabilization  of  a  mach  5,92  flat-plate  boundary  layer  using  local  sections  of 
porous  coatings 

Publications:  [3],  [9],  [13],  [15],  [20] 

Passive  eontrol  of  boundary  layer  transition  by  using  porous  eoating  to  stabilize  hypersonie 
boundary  layers  over  flat  plates  and  eones  has  been  sueeessfully  demonstrated  by  theoretieal 
analyses,  experiments,  and  numerieal  simulations  [14-18].  Fedorov  et  al.  [14]  performed  theoretieal 
analyses  on  the  seeond-mode  stability  of  a  hypersonie  boundary  layer  over  a  flat  plate  eovered  by 
an  ultrasonieally  absorptive  eoating  (UAC).  They  found  that  the  seeond  mode  growth  was 
massively  redueed,  beeause  the  porous  layer  absorbed  the  disturbanee  energy.  To  demonstrate  the 
stabilization  effeet  of  UAC  on  the  seeond  mode,  Rasheed  et  al.  [15]  experimentally  studied  the 
stability  of  a  Maeh  5  boundary  layer  on  a  sharp  5.06-deg  half-angle  eone  at  zero  angle  of  attaek. 
The  eone  had  a  smooth  surfaee  around  half  the  eone  eireumferenee  and  an  UAC  porous  surfaee  on 
the  other  half.  Their  experiments  indieated  that  the  porous  surfaee  was  highly  effeetive  in  stabilizing 
the  seeond  mode  and  delaying  transition,  when  the  pore  size  was  signifieantly  smaller  than  the 
disturbanee  wavelength.  Fedorov  et  al.  [16]  experimentally  and  theoretieally  studied  the  effeet  of  an 
UAC  on  hypersonie  boundary  layer  instabilities.  Their  experiments  were  performed  on  a  7-degree 
half-angle  sharp  eone  in  a  Maeh  6  wind  tunnel.  Half-surfaee  of  the  eone  was  solid  whereas  the  other 
half-surfaee  was  eovered  by  a  felt-metal  eoating.  Both  theoretieal  predietions  and  experimental 
measurements  on  two-  and  three-dimensional  disturbanees  showed  that  the  porous  eoating  strongly 
stabilized  the  seeond  mode  and  marginally  destabilized  the  first  mode.  Maslov  [17]  experimentally 
studied  the  stabilization  of  hypersonie  boundary  layer  by  miero  struetural  porous  eoating.  In  his 
experiments,  both  regular  porous  UAC  and  random  felt-metal  porous  UAC  are  used.  The  results 
eonflrmed  that  porous  eoating  strongly  stabilizes  the  seeond  mode  and  marginally  destabilizes  the 
first  mode.  Compared  with  regular  strueture  UAC,  felt-metal  UAC  is  mueh  stronger  in  first  mode 
destabilization,  with  the  peak  amplitude  of  the  first  mode  inereasing  around  70%.  Egorov  et  al.  [18] 
studied  the  effeet  of  porous  eoating  on  stability  and  reeeptivity  of  a  Maeh  6  flat-plate  boundary 
layer  by  two-dimensional  numerieal  simulation  using  a  seeond-order  TVD  seheme.  They  found  that 
a  porous  eoating  of  regular  porosity  effeetively  diminishes  the  seeond  mode  growth  rate,  while 
weakly  affeeting  aeoustie  waves. 


10 


In  previous  studies,  porous  eoating  eovers  either  the  entire  flat  plate  or  the  surfaee  around  half  the 
eone  eireumferenee.  The  effeet  of  porous  eoating  loeations  on  boundary  layer  stabilizations  has  not 
been  eonsidered.  Furthermore,  the  destabilization  effeet  of  porous  eoating  on  the  Maek  first  mode 
has  not  been  studied  in  detail.  However,  how  to  loeate  porous  eoatings  may  affeet  the  effleieney  of 
the  boundary  layer  stabilization.  The  first  mode  at  a  given  frequeney  is  spatially  developed  into  the 
seeond  mode  at  the  same  frequeney  when  propagating  downstream.  Therefore,  destabilization  of  the 
first  mode  deereases  the  effleieney  of  the  seeond  mode  stabilization.  In  eurrent  researeh,  the 
stabilization  of  a  Maeh  5.92  boundary  layer  over  a  flat  plate  using  loeal  seetions  of  porous  eoating 
is  studied  by  a  eombination  of  direet  numerieal  simulation  (DNS)  and  linear  stability  theory  (LST), 
foeusing  on  the  effeet  of  porous  eoating  loeations  and  the  destabilization  of  the  Maek  first  mode. 
The  stability  simulations  eonsist  of  two  steps:  1 .  periodie  disturbanees  eorresponding  to  mode  S  or 
mode  F  are  superimposed  on  steady  base  flow  at  a  eross-seetion  of  the  boundary  layer  to  show 
spatial  development  of  the  wave;  2.  loeal  seetions  of  porous  eoating  are  used  downstream  of  the 
superimposed  wave  to  investigate  its  effeet  on  boundary  layer  instability.  The  model  of  the 
felt-metal  porous  eoating  is  the  same  as  that  of  Fedorov  et  al.  [16].  At  step  2,  a  series  of  numerieal 
simulations  are  earned  out  by  putting  porous  eoatings  both  upstream  and  downstream  of  the 
synehronization  point.  Numerieal  simulation  results  are  interpreted  by  eomparing  with  the 
theoretieal  analyses  of  LST. 


Fig.  5  Spatial  developments  of  superimposed  mode  S  and  mode  F:  pressure  amplitude  and  growth  rates 

Fig.  6  shows  the  spatial  developments  of  superimposed  mode  S  and  mode  F.  Periodie  disturbanees 
eorresponding  to  mode  S  or  mode  F  at  a  frequeney  of  100  kHz  is  superimposed  on  steady  base  flow 
at  a  eross-seetion  of  the  boundary  layer  at  x*  =  69.00  mm  (R  =  987.24).  In  this  figure,  the  spatial 

o 

development  of  superimposed  mode  F  with  an  amplitude  parameter  being  10'  is  also  plotted.  It  is 
elearly  shown  that  the  spatial  developments  of  mode  F  with  different  amplitudes  have  the  same 
profde.  One  hundred  times  differenee  in  superimposed  disturbanee  amplitude  leads  to  100  times 
amplitude  differenee  in  spatial  development,  whieh  illustrates  the  linear  properties  of  disturbanees. 
After  mode  S  is  superimposed  at  x*  =  69.00  mm,  it  grows  dramatieally  with  the  peak  of  amplitude 
being  around  x*  =  0.8  m.  The  seenario  of  mode  F  is  quite  different  from  that  of  mode  S.  After  the 
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transient  process,  amplitude  of  mode  F  is  decreasing  because  mode  F  is  a  stable  wave.  The  lowest 
amplitude  of  mode  F  is  observed  at  around  x*  =  0.24  m.  Then  the  wave  grows  quite  similar  to  mode 
S,  also  with  the  peak  of  amplitude  being  around  x*  =  0.8  m. 

Fig.  7  also  compares  the  growth  rate  calculated  from  stability  simulations  with  that  obtained  from 
LST  for  superimposed  mode  S  and  mode  F.  For  mode  S,  the  growth  rate  of  stability  simulation  has 
a  good  agreement  with  that  of  LST  in  the  region  from  ro  =  0.1 1  to  ro  =  0.13.  When  ro  is  larger  than 
0.13,  the  growth  rate  of  stability  simulation  is  larger  than  that  of  LST,  i.e.,  mode  S  obtained  by 
stability  simulation  becomes  more  stable  than  that  predicted  by  LST.  For  mode  F,  the  figure  shows 
that  the  growth  rate  of  stability  simulation  initially  has  a  good  agreement  with  that  of  LST.  However, 
it  approaches  that  of  mode  S  near  the  synchronization  point,  which  indicates  that  mode  F  changes  to 
mode  S  near  the  synchronization  point.  After  the  change,  the  growth  rate  of  stability  simulation  has 
a  good  agreement  with  that  of  mode  S.  The  discrepancy  between  the  growth  rates  of  mode  S 
calculated  from  stability  simulation  and  LST,  when  ro  is  larger  than  0.13,  is  mainly  caused  by  the 
nonparallel  flow  effect. 


Fig.  8  Amplitude  distributions  of  pressure  perturbation  along  the  flat  plate. 

Fig.  9  shows  amplitude  distributions  of  pressure  perturbation  along  the  flat  plate  for  the  six  cases  of 
porous  coating,  together  with  the  spatial  development  of  mode  S.  It  is  noticed  that  pressure 
perturbations  of  all  the  six  cases  with  porous  coating  have  a  similar  profile  as  that  of  purely  mode  S 
without  porous  coating,  both  having  the  peak  of  amplitude  at  around  x*  =  0.8  m.  However,  porous 
coating  does  affect  the  amplitude  of  mode  S.  For  cases  1  to  3,  porous  coatings  are  upstream  of  the 
synchronization  point.  Amplitude  of  pressure  perturbation  increases  with  the  number  of  porous 
coating  increasing  from  case  1  to  case  3.  For  cases  4  to  6,  porous  coatings  except  the  first  three  are 
downstream  of  the  synchronization  point.  Amplitude  of  pressure  perturbation  decreases  with  the 
number  of  porous  coating  increasing  from  case  4  to  case  6.  The  results  show  that  porous  coating 
destabilizes  the  Mack  first  mode  whereas  it  stabilizes  the  Mack  second  mode.  All  the  six  cases  of 
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stability  simulations  show  larger  peak  amplitude  of  mode  S,  which  means  destabilization  of  the 
Mack  first  mode  is  quite  significant.  In  order  to  control  boundary  layer  transition,  it  is  very 
important  to  put  porous  coating  downstream  of  the  synchronization  point,  where  mode  S 
corresponds  to  the  Mack  second  mode. 


3.4  Collaborated  with  Prof,  Tumin  on  receptivity  of  the  Mach  5,92  flow  over  a  flat  plate  to 
wall  blowing-suction 

Publications:  [4],  [10],  [19] 

Direct  numerical  simulations  of  receptivity  in  a  boundary  layer  over  a  flat  plate  are  carried  out  with 
2D  perturbations  introduced  into  the  flow  by  periodic-in-time  blowing-suction  through  a  slot  [5]. 
The  perturbation  flow  field  was  decomposed  into  normal  modes  with  the  help  of  the  multimode 
decomposition  technique  based  on  the  spatial  biorthogonal  eigenfunction  system.  The 
decomposition  allows  filtering  out  the  stable  and  unstable  modes  hidden  behind  perturbations 
having  another  physical  nature.  Method  of  multiple  scales  (MMS)  is  used  to  investigate  the 
nonparallel  flow  effects. 


Fig. 7  Comparisons  of  growth  rates  and  pressure  perturbation  ampiitudes. 

Fig.  7  compares  the  growth  rates  and  pressure  perturbation  amplitude  obtained  from  numerical 
simulations  with  those  predicted  by  LST  and  MMS.  Both  figures  show  that  the  predictions  of  MMS 
have  a  good  agreement  with  the  numerical  results.  Nonparallel  flow  effects  destabilize  the  first 
mode  whereas  stabilize  the  second  mode.  Also  it  is  caused  by  the  nonparallel  flow  effects  that,  in 
the  numerical  simulation  results,  the  growth  rate  of  the  second  mode  is  always  smaller  than  that 
predicted  by  the  linear  stability  theory. 

3.5  Studied  the  nose-bluntness  and  entropy-layer  effects  on  hypersonic  boundary  layer  in  the 
context  of  transient  growth  with  or  without  surface  roughness. 
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Publications:  [16],  [24] 


For  hypersonic  boundary  layers  over  blunt  cones,  the  transition  reversal  phenomenon  refers  to  the 
experimental  observation  that  the  laminar-turbulent  transition  location  moves  upstream  with 
increasing  nose  radius  when  it  is  larger  than  a  certain  critical  value.  Currently,  there  is  no 
satisfactory  explanation  for  transition  reversal  phenomenon.  All  previous  linear  stability  theory 
(LST)  and  numerical  simulation  studies  on  the  nose  bluntness  effects  have  not  found  transition 
reversal.  The  previous  theoretical  and  computational  studies,  however,  have  not  been  done  on  the 
actual  experimental  flow  conditions  which  had  shown  transition  reversal. 

Surface  roughness  has  been  shown  to  have  profound  effects  on  boundary  layer  transition.  Recent 
research  has  shown  that  one  possible  explanation  to  bypass  transition  is  the  transient  growth  theory. 
However,  there  has  not  been  any  direct  numerical  simulation  study  on  transient  growth  in 
hypersonic  boundary  layers.  Furthermore,  it  is  not  known  how  the  optimal  disturbances  computed 
by  the  transient  growth  theory  are  generated  by  surface  roughness.  This  paper  presents  some  initial 
results  in  our  on-going  numerical  simulation  study  on  the  receptivity  of  the  hypersonic  boundary 
layers  to  stationary  surface  roughness  without  any  freestream  forcing  waves.  The  flow  conditions 
are  those  of  Stetson’s  1984  wind-tunnel  experiments  for  Mach  7.99  flow  over  a  7  degree  half-angle 
blunt  cone  [19].  The  simulation  results  show  that  the  roughness  element  generates  streamwise 
vortices  inside  the  boundary  layer.  In  the  computational  domain  studied  in  current  project,  the 
induced  perturbation  mainly  decays  after  the  roughness.  No  strong  transient  growth  appears 
immediately  behind  the  roughness  element. 


Fig. 8  Distributions  of  u  perturbation  magnitudes  for  aii  14  modes  aiong  a  streamwise  grid  iine  iocated  near  the 

edge  of  the  boundary  iayer. 

Fig.  8  shows  the  distribution  of  u  perturbation  magnitudes  for  all  14  modes  along  a  streamwise  grid 
line  located  near  the  edge  of  the  boundary  layer.  The  perturbations  rise  rapidly  for  flow  over  the 
roughness.  Downstream  of  the  roughness,  the  perturbations  drop  initially  for  all  n.  On  the  other 
hand,  the  perturbations  for  1  n  =  start  to  increase  after  the  initial  drop,  while  the  perturbations  for 
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other  modes  decrease.  No  obvious  strong  transient  growth  can  be  observed  in  the  current  simulation 
within  the  limited  length  of  the  flow  field  behind  the  roughness  element. 
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Fig. 9  Cross-sectional  view  of  the  contours  of  u  perturbations  for  two  of  the  14  modes. 


Fig.  9  shows  the  cross-sectional  view  of  the  contours  of  temperature  and  horizontal  velocity 
perturbations  for  two  of  the  14  modes  (n  =  1  and  8).  The  features  for  all  other  modes  are  similar  to 
those  in  the  figures.  For  the  cases  of  n  =  1  and  n  =  S  ,  Fig.  1 1  shows  that  the  flow  decelerates  behind 
the  peak  of  the  roughness  heights,  and  accelerates  behind  the  valleys  of  the  roughness  element,  as  a 
results  of  the  streamwise  vortices  generated  by  the  surface  roughness.  The  results  indicate  the 
generation  of  the  streamwise  vortices  by  the  surface  roughness. 


4  SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 


In  the  report  period,  we  have  mainly  focused  on  the  DNS  studies  of  receptivity  of  supersonic  and 
hypersonic  boundary  layer  flows.  A  new  high-order  cut-cell  method  has  been  developed  for  the 
numerical  simulation  of  hypersonic  boundary  layer  transition  with  finite  height  surface  roughness. 
The  method  has  been  applied  to  the  numerical  simulations  of  two-dimensional  hypersonic  flows 
over  a  fiat  plate.  Furthermore,  the  stabilization  effect  of  the  surface  porous  coating  over  a  flat  plate 
is  extensively  studied  by  series  of  numerical  simulations.  We  also  collaborate  with  Prof.  Tumin  in 
the  University  of  Arizona  to  compare  numerical  and  theoretical  results  on  receptivity  of  a  Mach 
5.92  flow  over  a  flat  plate  to  wall  blowing-suction,  and  to  deeply  analyze  the  nonparallel  flow 
effects.  The  most  signiffcant  research  accomplishments  are 

1.  We  developed  a  new  high-order  cut-cell  method  for  computing  boundary  layer  transition  with 
ffnite  surface  roughness.  Due  to  the  difficulty  in  grid  generation  for  numerical  simulation  of 
high-speed  flow  with  arbitrary  surface  roughness  elements,  it  is  advantageous  to  use  a  fixed-grid 
cut-cell  method  to  compute  such  flow.  Our  new  method  is  mainly  based  on  finite  difference 
method,  and  adopts  Gibou  and  Fedkiw’s  shifting- interpolation-stencil  strategy,  i.e.,  the  “drop 
point”  approach,  to  avoid  the  “smell  cell”  problem.  The  new  cut-cell  method  is  uniformly 
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high-order  for  all  grids  in  the  eomputational  domain.  The  eomputational  aeeuraey  of  our  third 
and  fourth  order  aeeurate  eut-eell  methods  is  tested  for  the  eomputations  of  Poisson  and 
hyperbolie  equation  with  Diriehlet  boundary  eonditions  in  irregular  domains.  We  then  test  a 
third-order  aeeurate  eut-eell  method  for  the  two-dimensional  Navier-Stokes  equations.  The  new 
seheme  is  used  to  simulate  steady  and  unsteady  hypersonie  boundary  layer  flows  over  a  flat 
plate  with  an  isolated  surfaee  roughness  element  (Fig.  1). 

2.  We  have  studied  the  stabilization  of  a  Maeh  5.92  over  a  flat  plate  to  two-dimensional  wall 
blowing-suetion  perturbation  with  finite  roughness  loeated  in  different  loeations.  The  numerieal 
results  show  that  the  development  of  two-dimensional  wall  perturbations  eventually  results  in 
the  same  type  of  instability  wave  (mode  S)  in  the  boundary  layer  with  roughness.  The  roughness 
loeation  plays  an  important  role  in  the  developments  of  mode  S  by  the  blowing-suetion  aetuator. 
Mode  S  is  destabilized  only  when  the  roughness  element  is  loeated  elose  upstream  of  the 
synehronization  point.  On  the  other  hand,  when  the  roughness  element  in  is  downstream  of  the 
synehronization  point,  mode  S  is  stabilized.  This  happens  even  when  the  roughness  is  still 
within  the  unstable  region  of  mode  S.  Therefore,  the  synehronization  point  and  finite  roughness 
loeation  is  eritieal  to  the  reeeptivity  proeess.  The  relationship  between  the  loeation  of  the 
roughness  and  the  synehronization  point  suggests  that,  in  order  to  eontrol  or  delay  the 
laminar-turbulent  transition  more  effieiently,  the  roughness  element  should  be  plaeed  upstream 
of  the  synehronization  point  between  mode  S  and  mode  F. 

5  DEVELOPMENT  OF  A  NEW  HIGH-ORDER  CUT-CELL  METHOD 


One  major  aeeomplishment  of  the  eurrent  researeh  projeet  is  the  development  of  a  new  high-order 
eut-eell  method.  Our  new  method  is  mainly  based  on  finite  differenee  method,  and  adopts  Gibou 
and  Fedkiw’s  shifting-interpolation-steneil  strategy,  i.e.,  the  “drop  point”  approaeh,  to  avoid  the 
“smell  eell”  problem.  The  new  eut-eell  method  is  uniformly  high-order  for  all  grids  in  the 
eomputational  domain.  The  eomputational  aeeuraey  of  our  third  and  fourth  order  aeeurate  eut-eell 
methods  is  tested  for  the  eomputations  of  Poisson  and  hyperbolie  equation  with  Diriehlet  boundary 
eonditions  in  irregular  domains.  We  then  test  a  third-order  aeeurate  eut-eell  method  for  the 
two-dimensional  Navier-Stokes  equations.  The  new  seheme  is  used  to  simulate  steady  and  unsteady 
hypersonie  boundary  layer  flows  over  a  fiat  plate  with  an  isolated  surfaee  roughness  element.  The 
new  high-order  eut-eell  method  has  been  tested  and  published  in  Publications:  [7],  [18]. 

5,1  Introduction 

The  understanding  of  physieal  meehanisms  of  roughness  indueed  boundary-layer  transition  is 
eritieal  to  the  development  of  hypersonie  vehieles  [20].  Transition  ean  have  a  first-order  impaet  on 
the  lift  and  drag,  stability  and  eontrol,  and  heat  transfer  properties  of  the  vehieles  [21].  For  example, 
roughness  indueed  transition  is  an  important  eonsideration  in  the  design  of  thermal  proteetion 
systems  (TPS)  of  hypersonie  vehieles  [22,  23].  Figure  10  shows  an  example  surfaee  roughness  on 
test  models  for  hypersonie  boundary  layer  transition.  For  a  reentry  vehiele  entering  earth’s 
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atmosphere,  it  initially  experienees  a  heating  environment  assoeiated  with  a  laminar  boundary  layer. 
As  the  vehiele  attitude  deereases,  the  vehiele  surfaee  beeomes  rougher  and  the  boundary  layer 
beeomes  turbulent.  The  transition  from  a  laminar  boundary  layer  to  a  turbulent  one  leads  to  the 
inerease  of  surfaee  heating  rates  by  a  faetor  of  five  or  more.  Thus  the  ability  to  understand  and 
prediet  the  physies  of  roughness  indueed  transition  plays  an  essential  role  in  the  design  of  TPS  for 
reentry  vehieles.  Currently,  surfaee  roughness,  espeeially  arbitrary  roughness  indueed 
laminar-turbulent  transition  in  hypersonie  boundary  layers  is  still  poorly  understood  due  to  the 
limitation  in  experimental  faeilities  and  numerieal  methods  [24]. 


Fig.  10.  An  example  surface  roughness  on  test  models  for  hypersonic  boundary  layer  transition:  1)  isolated 

roughness,  2)  distributed  roughness  [25]. 

Direct  numerical  simulation  has  become  an  effective  research  tool  for  studying  hypersonic 
boundary  layer  receptivity,  stability  and  transition  by  numerically  solving  the  time-dependent 
three-dimensional  Navier-Stokes  equations  for  the  temporally  or  spatially  evolving  instability 
waves.  It  is  necessary  to  use  high-order  numerical  methods  for  the  simulation  in  order  to  resolve  the 
wide  range  of  length  and  time  scales  of  the  complex  wave  fields  in  hypersonic  boundary  layers. 
Hence,  high-order  finite-difference  schemes  have  recently  received  much  attention  for  the  direct 
numerical  simulations  of  transitional  and  turbulent  flows  [26-31].  Zhong  [1]  presented  a  fifth-order 
upwind  finite  difference  shock  fitting  method  for  the  direct  numerical  simulation  of  hypersonic 
flows  with  a  strong  bow  shock  and  with  stiff  source  terms.  The  use  of  the  shock-fitting  method 
makes  it  possible  to  accurately  compute  the  shock-disturbance  interactions,  and  the  development  of 
instability  waves  in  the  boundary  layers.  The  fifth-order  shock-fitting  schemes  were  derived  on  a 
uniform  grid.  For  a  curvilinear  stretched  grid,  typically  used  in  simulations  of  viscous  flow  in  a 
boundary  layer,  the  physical  coordinates  with  a  non-uniform  grid  are  first  transformed  to  a  uniform 
grid  in  the  computational  coordinates.  The  high-order  schemes  are  subsequently  applied  to  the 
transformed  equations  in  a  uniform  grid.  The  fifth-order  shock-fitting  scheme  was  used  for 
numerical  studies  of  the  receptivity  of  two-dimensional  Mach  15  flows  over  a  blunt  leading  edge 
[32]. 

Most  of  the  DNS  methods,  however,  have  been  developed  for  smooth  surfaces  computed  by 
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body-fitted  grids.  Such  methods  may  not  be  able  to  compute  flow  with  surface  roughness  elements 
similar  to  those  shown  in  Fig.  10.  The  main  obstacle  is  the  difflculty  in  generate  smooth  body-fitted 
grids  around  the  surface  of  arbitrary  roughness.  One  approach  for  overcoming  the  difficulty  in  grid 
generation  is  to  use  a  Cartesian  grid  method,  which  is  easy  in  grid  generation.  Cartesian  grid 
methods  can  also  take  full  advantage  of  fast  computer  architectures  such  as  parallel  computers  and 
can  serve  as  a  very  flexible  method  for  simulating  flow  around  complex  geometries.  On  the  other 
hand,  a  Cartesian  grid  method,  the  grid  lines  are  not  aligned  with  the  body  surface.  As  a  result, 
special  treatment  is  needed  to  compute  the  flow  equations  in  local  regions  adjacent  to  the  solid 
surface. 

Various  Cartesian  grid  methods  have  been  developed  and  used  to  solve  problems  with  arbitrary 
geometry  [33-37].  Peskin  developed  the  immersed  boundary  method  (IBM)  [37]  for  the  simulation 
of  blood  flow  in  hearts.  The  method  is  based  on  Cartesian  grids  where  the  surface  of  solid 
membrane  is  represented  by  a  discrete  delta  function.  This  delta  function  is  added  into  the 
Navier-Stocks  equations  to  reflect  the  immersed  boundary  between  the  two  phase  flows.  The 
resulting  equations  are  discretized  by  a  standard  finite  difference  method  in  a  fixed  Cartesian  grid 
system.  Since  its  introduction,  the  IBM  methods  have  been  applied  to  many  different  fluid  flow 
problems,  including  flow  interaction  with  solid  surface.  For  example,  Marxen  et  al.  [38]  applied  the 
IBM  to  simulate  the  effects  of  a  localized  two-dimensional  roughness  element  on  the  disturbance 
amplification  in  a  hypersonic  boundary  layer.  However,  since  the  immersed  boundary  method  uses 
the  discrete  delta  function  approach,  it  leads  to  a  smeared  interface  with  a  thickness  in  an  order  of  a 
mesh  width.  The  immersed  boundary  method  is  locally  first-order  accurate  at  the  interface,  which 
may  not  be  accurate  enough  for  the  DNS  of  hypersonic  boundary  transition  problems. 

In  contrast  to  the  immersed  boundary  method  with  smeared  interfaces,  “sharp  interface”  Cartesian 
grid  methods,  which  maintain  second-order  accuracy  at  the  interfaces,  have  been  developed  [34,  39]. 
Udaykumar  et  al.  [34,  35]  used  the  finite-volume  methodology  to  solve  the  incompressible 
Navier-Stokes  equations  for  flow  interacting  with  moving  bodies.  The  method  is  based  on  a  fixed 
Cartesian  mesh  where  the  solid  boundaries  can  move  across  the  grid  lines.  The  flow  equations  are 
discretized  by  an  overall  second-order-accurate  finite-volume  technique.  The  interface  is 
represented  as  a  sharp  boundary  between  the  fluid  and  solid  phases.  A  one-sided  bilinear 
interpolation  is  applied  to  calculate  the  accurate  flow  conditions  in  the  sharp  irregular  boundary. 
Johansen,  McCorquodale  and  Colella  [39,  40]  developed  a  similar  sharp-interface  Cartesian  grid 
method  for  solving  two-dimensional  Poisson  and  heat  equations  on  irregular  domains.  Their  method 
was  based  on  finite-volume  formulation,  with  imbedded  irregular  grid  cells  on  the  boundary.  The 
irregular  grid  cells  were  treated  by  conservative  differencing  of  second-order  accurate  fluxes  on 
each  cell  volume.  The  method  is  second  order  accurate  for  a  problem  with  irregular  boundary. 

Fedkiw  and  Gibou  [33,  41]  developed  a  ghost  fluid  method,  which  is  a  different  implementation  of 
the  sharp  interface  methods,  for  solving  flow  equations  in  arbitrary  domain.  The  basic  idea  was  to 
extrapolate  variables  on  one  side  of  the  boundary  into  a  number  of  ghost  cells  located  on  other  side 
of  the  solid  boundary.  The  method  was  uniformly  second  order.  Gibou  and  Fedkiw  [41]  improved 
accuracy  of  the  ghost  fluid  method  to  fourth-order  accuracy  for  both  the  Laplace  equation  and  the 
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heat  equation  with  Diriehlet  boundary  eonditions  on  irregular  domains.  Finite  difFerenee 
diseretization  was  used  to  solve  the  equations  in  an  irregular  domain.  The  sharp  interfaee  Cartesian 
grid  method  has  also  been  termed  eut-eell  method  in  handling  irregular  grid  eells  along  the 
boundary  [42]. 

For  multi-phase  flows  with  moving  interfaees,  an  immersed  interfaee  method  was  developed  and 
has  been  widely  used  [43-53].  For  an  immersed  interfaee  method,  a  Cartesian  grid  is  often  used 
where  a  sharp  interfaee  moves  freely  aeross  the  fixed  grid  lines.  This  method  ean  aehieve  a  seeond 
or  higher  order  global  aeeuraey  by  ineorporating  jump  eonditions  into  the  finite  differenee  formulas 
for  variables  and  their  derivatives  at  the  interfaee.  This  method  was  first  developed  by  Leveque  and 
Li  [12]  to  solve  elliptie  equations  with  diseontinuous  eoeffieients  and  singular  sourees.  The  original 
method  was  seeond  order  in  the  interior  of  the  domain  and  loeally  first  order  at  the  interfaee. 
Leveque  and  Li  [54]  subsequently  extended  the  method  to  the  Stokes  flow  with  elastie  boundaries 
or  surfaee  tension.  Wiegmann  and  Bube  [46,  47]  developed  an  explieit  jump  immersed  interfaee 
method  for  speeial  eases  where  the  explieit  jump  eonditions  of  variables  and  their  high-order 
derivatives  are  known.  This  method  ean  only  aehieve  arbitrary  high-order  aeeuraey  if  the 
eorresponding  high-order  derivatives  of  jump  eonditions  ean  be  analytieally  derived.  Zhong  [53] 
developed  a  new  high-order  immersed  interfaee  method  whieh  does  not  require  seeond  and  higher 
derivatives  of  jump  eonditions.  The  main  idea  was  to  use  a  wider  grid  steneil  aeross  the  interfaee 
instead  of  taking  higher  derivatives  of  jump  eonditions  to  aehieve  high-order  aeeuraey  for  finite 
differenee  formulas  at  the  interfaee. 

For  sharp-interfaee  Cartesian  grid  methods,  a  “small  eell  problem”  [55]  of  numerieal  instability 
would  arise  when  finite-differenee  or  finite-volume  methods  are  applied  to  relatively  small-sized 
irregular  grid  eells  ereated  by  a  sharp-interfaee  Cartesian  grid  method.  For  a  time-dependent 
problem,  the  small  eell  problem  will  signifieantly  restriet  the  size  of  times  steps  in  temporal 
integration  methods.  Many  methods  have  been  proposed  to  resolve  the  small  eell  problem.  Berger 
and  Leveque  [55]  used  a  rotating  box  method.  Johansen  and  Colella  [39]  used  a  flux-redistribution 
proeedure.  Quirk  [56]  and  Udaykumar  [35]  employed  a  eell  merging  method  to  avoid  small  eells  in 
order  to  maintain  numerieal  stability. 

Gibou  and  Fedkiw  [41]  presented  a  stable  fourth-order  finite  differenee  method  for  solving  the 
Laplaee  equation  on  an  irregular  domain.  They  showed  that  a  good  rule  of  thumb  for  removing 
small  eell  restrietion  is  that  the  interpolation  should  be  shift  to  be  eentered  one  grid  points  left  or 
right  if  the  interpolation  steneil  involves  point  eonsidered  to  be  relatively  elose  to  the  boundary. 

Due  to  the  diffieulty  in  grid  generation  for  numerieal  simulation  of  high-speed  flow  with  arbitrary 
surfaee  roughness  elements,  it  is  advantageous  to  use  a  fixed-grid  eut-eell  method  to  eompute  sueh 
flow.  The  grid  ean  be  a  smooth  eurvilinear  body-fitted  one  along  a  baseline  smooth  surfaee  without 
the  roughness.  The  aetual  surfaee  with  roughness  will  eut  aeross  the  grid  lines,  whieh  ereate 
irregular  grid  points.  Most  of  the  previous  Cartesian-grid  eut-eell  (or  sharp-interfaee)  methods  with 
different  boundary  treatments,  however,  are  only  first  or  seeond  order  aeeurate  at  the  interfaee.  The 
relative  low  aeeuraey  at  the  interfaee  may  not  be  suffieient  for  numerieally  simulating 
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laminar-turbulent  boundary-layer  transition  involving  surfaee  roughness.  For  a  simulation  study  of 
hypersonie  boundary  layer  transition  shown  in  Fig.  1,  it  is  neeessary  to  eonsider  the  effeets  of  both  a 
bow  shoek  and  an  arbitrary  roughness  element  on  the  lower  surfaee.  More  details  ean  be  found  in 
[1].  Numerieal  methods  for  sueh  studies  need  to  be  high-order  aeeurate  both  in  resolving  the 
eomplex  shoek/disturbanee  interaetion  and  in  resolving  a  wide  range  of  flow  time  and  length  seales. 
In  addition,  high-order  aeeuraey  is  required  for  the  loeal  finite  differenee  sehemes  for  loeal  region 
adjaeent  to  the  surfaee  of  the  roughness  element.  In  this  paper,  we  present  a  new  high-order  eut-eell 
method  for  numerieal  simulation  of  roughness  indueed  transition  in  viseous  hypersonie 
boundary-layer  flows.  This  method  is  mainly  based  on  finite  differenee  method,  and  is  different  in 
implementation  from  all  the  methods  diseussed  above.  To  overeome  the  small  eell  problem,  we 
adopt  Gibou  and  Fedkiw’s  shifting-interpolation-steneil  strategy,  i.e.,  the  “drop  point”  approaeh. 
The  new  eut-eell  method  is  uniformly  high-order  for  both  regular  and  irregular  grid  points.  The 
eomputational  aeeuraey  of  our  third  and  fourth  order  aeeurate  eut-eell  methods  is  tested  for  the 
eomputations  of  Poisson  and  hyperbolie  equation  with  Diriehlet  boundary  eonditions  in  irregular 
domains.  We  then  test  a  third-order  aeeurate  eut-eell  method  for  the  two-dimensional  Navier-Stokes 
equations.  The  new  seheme  is  used  to  simulate  steady  and  unsteady  hypersonie  boundary  layer 
flows  over  a  flat  plate  with  an  isolated  surfaee  roughness  element  (Fig.  1).  The  bow  shoek  generated 
from  the  leading  edge  of  the  flat  plate  is  treated  as  a  eomputational  boundary  and  diseretized  based 
on  Zhong’s  [1]  high-order  finite-differenee  shoek-fitting  method.  The  unsteady  flow  over  flat  plate 
is  exeited  by  a  narrow  blowing  and  suetion  slot  eomposed  of  15  different  frequeneies  mounted  near 
the  leading  edge.  Several  multi-dimensional  eases  are  tested  to  demonstrate  the  order  of  aeeuraey 
for  these  numerieal  methods. 

5,2  High-Order  Cut-Cell  Method 

Our  high-order  eut-eell  method  is  presented  in  this  seetion  for  the  numerieal  simulation  of 
hypersonie  boundary-layer  transition  with  arbitrary  surfaee  roughness.  The  traditional  numerieal 
methods  of  ehoiee  for  the  DNS  of  transitional  and  turbulent  flows  have  been  speetral  methods 
beeause  of  their  high  aeeuraeies  [57-59].  But  the  applieations  of  speetral  methods  have  been  limited 
to  flows  in  simple  domains.  Finite-differenee  methods  have  reeently  reeeived  mueh  attention  for  the 
DNS  of  transitional  and  turbulent  flows,  espeeially  eompressible  flows,  [29,  60-63]  beeause  they 
ean  be  easily  applied  to  eomplex  geometries.  High-order  sehemes  are  required  beeause  traditional 
seeond-order  sehemes  do  not  provide  adequate  aeeuraey  level  for  the  direet  numerieal  simulation. 
Most  high-order  finite-differenee  methods  used  in  direet  numerieal  simulation  are  eentral  differenee 
sehemes  [29,  63]  whieh  introduee  only  phase  errors  but  no  dissipative  errors  in  numerieal  solutions. 
The  shortfall  of  eentral  sehemes  is  that  they  are  often  not  robust  enough  in  eonveetion  dominated 
hypersonie  flow  simulations.  On  the  other  hand,  Rai  et  al.  [60]  show  that  upwind-bias  sehemes  are 
very  robust  in  hypersonie  flow  simulation  even  when  they  are  made  high-order  aeeurate.  They  use  a 
spatially  fifth-order  upwind  finite-differenee  seheme  in  an  upwind-bias  steneil  to  eompute  the 
Navier-Stokes  equations.  The  numerieal  dissipation  in  the  upwind-bias  sehemes  is  enough  to 
eontrol  the  aliasing  errors  and  eould  maintain  the  overall  stability  of  the  method.  [64] 

5,2.1  Computational  Grids  and  Classification  of  Grid  Points 
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A  schematic  of  a  computational  domain  and  a  cut-cell  grid  in  roughness  induced  hypersonic 
boundary  layer  transition  is  shown  in  Fig.  11.  This  figure  shows  a  typical  hypersonic  flow  over  a 
blunt  body,  where  a  bow  shock  is  created  by  the  supersonic  freestream.  In  this  paper,  a  high-order 
shock-fitting  method  is  used  to  track  the  movement  of  the  bow  shock  which  is  treated  as  the  upper 
boundary  of  the  computational  domain.  The  computational  grid  for  a  shock  fitting  formulation  is 
bounded  between  the  bow  shock  above  and  the  blunt  body  below. 

The  cut-cell  grid  is  a  smooth  curvilinear  grid  fitted  to  the  baseline  body  shape  without  the  roughness. 
As  a  result,  the  roughness  surface  cuts  across  the  grid  lines.  The  roughness  surface,  F ,  is 
represented  by  surface  equation  in  the  following  form, 

r:/(x,y,z)  =  0  (1) 


For  a  problem  concerning  practical  arbitrary  roughness  in  hypersonic  vehicle  surface,  it  is  likely 
that  there  is  not  an  analytical  equation  applicable  to  represent  the  shape  of  the  roughness  element.  In 

this  case,  a  set  of  n  discrete  coordinate  points  along  the 

physical  roughness  surface  are  used  to  represent  the  surface.  With  these  roughness  data,  we  can 
reconstruct  the  roughness  surface  for  simulation  by  using  a  high-order  piecewise  polynomial 
interpolation. 


Fig.  11 .  Physical  and  computational  domain  and  a  cut-cell  grid  of  hypersonic  flow  over  a  blunt  body  with 
surface  roughness:  1)  physical  grid,  2)  computational  grid  with  a  transformed  roughness. 


Both  the  governing  Equation  and  the  roughness  equation  (1)  in  the  physical  domain  are  transformed 
into  a  Cartesian  computational  domain  bounded  by  bow  shock  and  fiat  plate  boundary.  Under  the 
computational  coordinate  system,  the  body  fitted  grids  are  represented  by  a  curvilinear 
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three-dimensional  eoordinates  {^,7],g,T)  along  the  grid  lines.  The  unsteady  movement  of  the  bow 


shoek  is  treated  as  the  eomputational  upper  boundary  loeated  at  ;/  =  ,  whieh  is  time  dependent. 

The  other  grid  lines  ^  =  const  and  g  =  const  remains  stationary  during  eomputations.  The 
eoordinate  transformation  is  defined  by: 

'^  =  ^(x,y,z) 
ri  =  ri{x,y,z,t) 

< 

g  =  g(x,y,z) 

T  =  t 


X  =  x(4,7J,g,T) 
y  =  y(i,fl,g,T) 
z  =  z(<^,/l,g,T) 
t  =  T 


(2) 


where  (x,y,  z,  t)  are  the  physieal  eoordinates  defined  under  Cartesian  eoordinate  system. 

Substituting  Eq.  (2)  into  the  governing  N-S  equation,  we  obtain  a  system  of  transformed  governing 
equations  in  the  eomputational  domain  {^,ri,g,T)  as 


\dU  BE'  BE'  BG'  BE'  BE'  BG'  ,,5(1//) 

■  + - + - + +  — ^  +  +  — —  +  U— -  =  0 


J  Bt  B^  Btj  B(^  B^  Btj  B(^  Bt 


The  transformed  fluxes  of  the  equation  above  are: 


J 


J 


^  X  ^2^ y  ^2^ z 

J 


J 


e:  = 


J 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
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G. 


J 


(9) 


are  transformation  metrics,  and  J  is  Jacobean  matrix  of  coordinate 

transformation  defined  by 


S{^,ri,g) 

d{x,y,z) 


(10) 


In  addition  to  the  transformation  of  the  governing  equations,  the  equation  for  the  surface  roughness 
is  also  transformed  into  the  computational  domain  and  can  be  represented  as 

y,  z),7]{x,  y,  z,  t),  g{x,  y,  z))  =  0  (11) 


The  grids  transformation  metrics  and  the  Jacobean  matrix  J  can  be  obtained  either  by  analytical 
formulas  of  the  coordinate  transformation  or  by  numerical  approximation. 

With  the  coordinate  transformation,  a  set  of  uniformly  distributed  Cartesian  grids  can  be  generated 
in  the  computational  domain  where  the  grid  distribution  in  the  physical  domain  is  not  uniformly 
distributed.  Because  smooth  body-fitted  grids  are  generated  in  the  regular  computational  domain 
without  the  roughness  as  shown  in  Fig.  11  (b),  some  of  the  Cartesian  grid  cells  may  be  cut  by  the 
roughness  boundary,  which  leads  to  irregular  Cartesian  grid  cells. 
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Fig.  12.  Classification  of  four  types  of  grids  for  the  high-order  cut-cell  method  in  computational  domain  (1^,7) , 
where  V  represents  regular  point,  ®  dropped  point,  O  boundary  point,  and  •  irregular  point. 
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In  the  current  high-order  cut-cell  finite-difference  method,  four  different  types  of  grid  points  are 
defined  according  to  their  relative  locations  with  respective  to  the  roughness  surface.  Different 
numerical  algorithms  are  implemented  for  different  types  of  points.  The  four  types  of  points,  as 
shown  in  Error!  Reference  source  not  found.  12,  are  defined  as: 

•  Boundary  points.  The  points  created  by  the  intersection  of  roughness  interface  F  and  grid 

lines  are  boundary  points.  They  are  not  part  of  the  original  Cartesian  grid,  and  they  are  used 
in  finite  difference  formulas  for  grid  points  adjacent  to  the  roughness  surface. 

•  Irregular  points:  For  those  points  located  close  enough  to  the  roughness  boundary  F  that 

their  finite  difference  stencils  contain  a  boundary  point,  they  are  defined  as  irregular  points. 
In  the  finite  difference  approximation  involving  an  irregular  point,  local  grid  spacing 
becomes  non-uniform  because  of  the  inclusion  of  a  boundary  point  in  the  stencil. 

•  Dropped  points.  If  a  grid  point  is  adjacent  to  a  boundary  point  along  a  grid  line  with  a  distance 

smaller  than  a  pre-specified  critical  ratio  0  in  the  ^  ox  rj  direction,  it  is  defined  as  a 

dropped  point  along  that  direction.  A  dropped  point  is  removed  from  the  grid  stencil  in  the 
formulation  of  a  local  finite  difference  approximation  in  the  corresponding  direction.  It 
should  be  noted  that  a  grid  point  can  be  a  drop  point  in  one  direction,  but  a  regular  point  in 
another  direction.  Therefore,  dropped  points  are  removed  in  the  “dropped  direction”  only. 
The  same  points  may  be  included  in  finite  difference  stencils  in  the  other  directions  if  these 
points  are  not  defined  as  “dropped”.  The  value  of  non-dimensional  critical  ratio  0  is  an 
adjustable  parameter.  In  the  cases  of  third  and  fourth  order  methods,  0  is  selected  to  be 
0.5  and  1.0  respectively.  All  of  the  points  on  the  solid  side  of  the  computational  domain 
are  defined  as  dropped  points  as  well.  They  do  not  participate  in  any  numerical  calculation. 

•  Regular  points:  All  other  grid  points  produced  by  the  intersection  of  grid  lines  themselves  are 

defined  as  regular  points.  Since  they  are  relatively  far  away  from  boundary  points,  a  standard 
finite  difference  approach  in  a  uniform  grid  can  be  applied. 


5.2.2.  Finite-Difference  Algorithms  for  Different  Types  of  Grid  Points 

The  derivatives  of  the  flux  terms  in  Eq.  (3)  are  discretized  by  different  methods  for  the  four 
different  types  of  grid  points.  The  flux  terms  in  regular  points  are  computed  by  the  standard  upwind 
finite  difference  scheme  introduced  in  next  section.  To  calculate  the  flux  in  irregular  points,  a 
high-order  non-uniform-grid  finite  difference  method,  the  stencil  of  which  consists  of  regular, 
irregular  and  boundary  points.  The  dropped  points  shown  in  Error!  Reference  source  not 
found.  12  are  not  included  in  non-uniform  or  uniform  grid  stencils  for  computing  flux  terms  in  a 
“dropped”  direction.  The  main  purpose  of  defining  a  dropped  point  is  to  avoid  the  small  cell 
problem  introduced  in  Section  1.  The  removal  of  drop  points  from  finite  difference  formulas 
ensures  that  the  distance  between  two  adjacent  grid  points  in  a  finite-difference  stencil  is  large 
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enough  in  order  to  prevent  from  the  deterioration  of  eorresponding  high-order  method  for 
eomputing  the  hyperbolie  part  of  Eq.  (3).  Based  on  the  different  grid  elassifieations  in  the 
eomputational  domain,  the  details  of  the  different  diseretization  methods  for  the  four  types  of  points 
are  deseribed  as  follows. 

5,2.2. 1  Regular  points  and  shock-fitting  methods 

Sinee  regular  points  are  not  affeeted  by  the  presenee  of  roughness  boundary,  regular  high-order 
upwind  sehemes  of  Zhong  [1]  are  used  for  the  diseretization  of  the  governing  equations  for  these 
points.  In  addition,  the  bow  shoek  shown  in  Fig.  11  is  treated  by  a  high-order  shoek-fitting  approaeh. 
The  finite  diflferenee  methods  and  formulas  used  for  the  regular  points  and  the  bow  shoek  are  briefly 
deseribed  here.  More  details  ean  be  found  in  Zhong  [1]. 

5,2,2,2.  Irregular  points 

Contrary  to  the  standard  fifth-order  finite-diflferenee  sehemes  used  for  regular  points,  speeial 
treatment  is  needed  in  the  diseretization  of  the  governing  equations  for  the  irregular  points  beeause 
boundary  points  are  ineluded  into  the  finite  differenee  steneil.  Fig.  13  shows  a  sehematie  of  a  grid 

steneil  for  irregular  points  near  the  boundary.  In  this  figure,  grid  points  ^2  ^3  are  irregular 

points,  4  is  a  boundary  point,  while  the  grid  point  between  and  ^2  ^  dropped  point 

beeause  it  is  too  elose  the  boundary  point.  The  dropped  point  is  removed  from  the  steneil  for  the 
irregular  point  sehemes.  The  rest  of  the  grid  points  are  regular  points.  For  example,  for  a  loeal 

third-order  finite-differenee  approximation  at  the  irregular  point  ^2 » the  grid  steneil  eonsists  of  the 

following  five  grid  points:  ,  ^2  >  •••’  •  The  grid  spaeing  between  the  points  involved  are  not 

uniform  beeause 


(j  = 


(12) 


where  9  is  the  spaeing  between  and  ^2  >  ls.h  is  the  uniform  grid  spaeing  of  the  regular 
grid.  Therefore,  a  non-uniform-grid  finite  diflferenee  sehemes  are  needed  flor  irregular  points. 
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Fig.  13.  A  schematic  of  a  gricJ  stencil  for  an  irregular  point  with  p  —  ^  ,  ^  =  5  and  Qf  =  {<^2,(^3}  ,  9  is 

non-uniform  grid  spacing  after  removing  a  dropped  point,  A/;  is  the  normal  grid  spacing,  0  represents  the 
dropped  point,  O  represents  the  boundary  point,  and  •  represents  irregular  and  regular  points. 

It  is  assumed  that  p  is  the  loeal  order  of  aeeuraey  in  boundary  of  the  simulation.  In  order  to 

maintain  a{p  +  \)-th  order  global  aeeuraey  for  the  distretization  of  both  the  inviseid  and  viseous 

flux  terms,  all  loeal  non-uniform  sehemes  for  irregular  points  need  to  be  at  least  p-th  order 

aeeuraey.  In  our  eonstruetion  of  high-order  eut-eell  seheme,  if  a  grid  steneil  with  a  number  of  q 
grids  is  used  for  diseretizing  regular  points  near  the  boundary  in  eaeh  direetion  of  the  eomputational 
domain,  there  are  L^/2J  irregular  points  near  a  boundary  surfaee.  We  denote  the  eolleetion  of  all 

irregular  points  near  this  boundary  in  one  direetion  as  set  QJ ,  where  the  superseript  represents  the 
direetion  of  the  steneil  and  the  subseript  represents  the  loeal  order  of  sehemes  at  this  point. 

As  a  general  eonvention,  the  boundary  point  is  labeled  as  ,  the  irregular  points  are  labeled  as 
^^‘l^^’^tially  in  the  order  of  their  distanees  from  the  boundary  point.  A  speeial  ease 

with  p  =  'i  and  q  =  5  is  shown  in  Fig.  13.  There  may  or  may  not  be  a  drop  point  in  the  steneil  for 

an  irregular  point  depending  on  the  spaeing  between  the  boundary  point  and  its  eloset  regular  grid 
point.  As  diseussed  earlier  in  this  paper,  if  the  non-dimensionalized  grid  spaeing  between  the 
dropped  point  0  and  the  boundary  point  O  is  less  than  a  pre-deseribed  eritieal  ratio  0  ,  this 
dropped  point  is  removed  from  the  grid  steneil  in  eorresponding  direetion.  Otherwise,  there  is  no 
dropped  point  in  the  steneil. 

Finite  difference  formulas  for  viscous  flux  terms  in  irregular  points 
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The  discretization  of  viscous  term  for  an  irregular  point  involves  the  calculations  of  second 

order  derivatives  since  contains  gradient  terms,  such  as  Vu  and  VT.  Under  the  coordinate 
transformation  of  Eq.  (2),  Vu  and  VT  can  be  expressed  as  derivatives  with  respect  to  coordinates  in 
the  computational  domain  {^,ri,g,T) .  A  two-step  algorithm  is  employed  for  the  computations  of 

viscous  fluxes.  First,  the  gradient  of  velocity  and  temperature  Vu  andVT  is  computed  at  each 
grid  point.  Using  these  results,  the  values  of  flux  vector  F^  is  evaluated  at  all  grid  points.  Second, 

the  derivatives  of  F^  are  computed  by  a  finite  difference  scheme.  Since  the  computations  in  both 

steps  involve  only  calculations  of  first-order  derivatives,  the  same  high-order  finite-difference 
schemes  are  used  at  each  step.  In  order  to  do  this,  the  first  derivative  has  to  be  computed  at  the 
boundary  points.  For  boundary  points  with  regular  finite-difference  stencil,  the  first-derivative  can 
be  calculated  by  using  one-side  finite-difference  scheme,  e.g.  for  4*'’  order  simulation: 


(wJi 


dWj  -  6^3  +  +  M5 

Ah 


(12) 


For  boundary  points  with  irregular  finite-difference  stencil,  the  first-derivative  can  be  evaluated  by 
using  non-uniform  one-side  finite-difference  scheme.  The  derivation  details  can  follow  the 
following  steps. 

The  general  formulation  of  a  non-uniform-grid  finite-difference  scheme  for  computing  the  viscous 
terms  for  an  irregular  point  of  QJ  can  be  written  in  the  following  form: 


^dF'^ 


Ji 


Ah  k=\ 


C{Ahy-^  idF 


q\ 


5^ 


for  e  where  i  =  2,3,...,[^ /2j  + 1  (13) 


The  subscripts  of  coefficient  a.  ^  stand  for  the  k-th  coefficient  for  the  i  -  th  irregular  point  as 

defined  for  the  case  of  p-l>  and  ^  =  5  (Fig.  13).  The  coefficient  a.  ^  is  function  of  <7  ,  which  is 
defined  in  Eq.  (12),  and  C  is  a  constant. 

The  coefficients  of  non-uniform-grid  finite-difference  equation  (Eq.  (13))  for  every  i-th  irregular 

grid  point  in  set  can  be  derived  either  by  a  Taylor  series  expansion,  or  by  taking  a  derivative  of 

a  polynomial  interpolated  through  the  non-uniform  stencil.  We  use  polynomial  interpolation  in  this 
paper.  Specifically,  for  each  irregular  point,  the  Fagrange  interpolation  polynomial  can  be  written  as 
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(14) 


q  (  q  f  ^  , 

/=1  \^m=l,mW  rm  y 

where  is  a  polynomial  interpolating  through  the  grid  steneil.  Differentiating  Eq.  (14)  onee 

with  respeet  to  ,  we  have 


^  ?  f  q  P—P  \  1 

=  V  V  TT  h _ i _  J7 

^  Lu  11  P —P  P —P 

l=\  \  n=\,n^l\  m=\,m^Lm^nhl  hm  )hl  J 


Substituting  ^  into  the  equation  above  and  eomparing  terms  with  those  of  Eq.  (13),  we 
obtain  the  eoeffieients  a.  ^  as  follows 

«a=  i  n  ^  ^  for  ^  =  1,2..^  (16) 

n=\,n*kym=l,m*k,m*n'Tk  hm  )  ~ 

where  a.  ^  is  a  funetion  of  cr  defined  in  Eq.  (12).  In  a  shoek-fitting  ealeulation,  the  eomputational 

grid  ehanges  with  the  movement  of  the  bow  shoek.  As  a  result,  the  value  of  6  is  a  funetion  of  time 
beeause  of  the  shoek  movement. 

During  the  time  advaneement  in  the  simulation,  the  grid  metries  in  the  eomputational  domain 
bounded  by  the  moving  shoek  vary  from  one  time  step  to  the  next.  As  a  result,  the  a  values  are  not 

eonstant.  Thus  we  need  to  reealeulate  the  finite-differenee  eoeffieients  a,.  ^  for  all  irregular  points 

in  eaeh  time  step.  In  our  simulation,  the  formulas  for  the  eoeffieients  are  stored  in  the  eomputer 
memory,  and  their  values  are  eomputed  explieitly  when  the  finite-differenee  sehemes  are 
implemented  in  eaeh  direetion. 

Finite  difference  formulas  for  inviscid  flux  terms  in  irregular  points 

Similar  to  the  loeal  algorithms  for  the  viseous  terms  at  irregular  points,  a  non-uniform-grid 
high-order  upwind  seheme  is  used  to  diseretize  the  inviseid  fluxes,  and  F'A  For  every 
irregular  grid  point,  there  are  several  possible  grid  steneils  for  finite-differenee  approximation  of  the 
flux  derivatives  of  the  same  aeeuraey  order.  Different  ehoiees  of  steneils  for  these  boundary  elosure 
sehemes  lead  to  different  stability  eharaeteristies  for  the  overall  algorithm.  Sinee  F'^  and  F'^ 
have  either  all  positive  or  all  negative  eigenvalues,  loeal  grid  steneils  for  finite-differenee 
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approximation  of  the  flux  derivatives  are  ehosen  so  that  the  diseretization  for  the  fluxes  are  upwind 
biased,  while  maintaining  high-order  aeeuraey.  Therefore,  for  a  given  irregular  grid  point,  the 
steneil  for  F'^  may  be  different  from  the  steneil  for  F'F  In  addition,  a  grid  steneil  for  these 
fluxes  may  or  may  not  inelude  the  boundary  point. 

The  non-uniform-grid  finite-differenee  sehemes  for  the  inviseid  flux  terms  of  positive  and  negative 
eigenvalues  ean  be  written  as 


1  1 

—  T  F'^ 

Ah  k=i 


for  e  where  i  =  2,3,...,\_q  /  2j  +  \ 


(17) 
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for  e  where  i  =  2,3,...,\_q 1 2\  +  \ 
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The  steneil  eontains  a  total  of  q  grid  points  as  shown  in  Fig.  13.  The  upwind  sehemes  are 
represented  by  different  sets  of  eoeffieients  of  the  two  formulas  above.  Beeause  the  upwind  bias 
dF'^ 

steneil  used  for  -  does  not  inelude  the  boundary  point,  the  eoeffieients  is  not  a  funetion 

of  cr  defined  in  Eq.  (12).  On  the  other  hand,  is  a  funetion  of  a  beeause  the  upwind  steneil 
dF'^ 

for  -  ineludes  the  boundary  point.  The  eoeffieients  for  eaeh  upwind  seheme  above  ean  be 

ealeulated  following  the  same  general  formulas  given  by  Eqs.  (15)  and  (16). 

The  higher-order  non-uniform  finite-differenee  steneils  require  suffieiently  many  grid  points  near 
the  roughness  surfaee.  When  the  grids  are  insuffieient  for  loeal  high-order  finite-differenee  steneil, 
there  are  two  approaehes  to  overeome  the  problem.  The  first  approaeh  is  to  maintain  the  original 
number  of  grids  point  in  the  loeal  high-order  finite-differenee  steneil.  More  grid  points  ean  be 
elustered  or  generated  near  the  irregular  boundary  by  either  adjusting  the  eoordination 
transformation,  or  refining  the  grids  in  the  eomputational  domain.  The  global  grids  refinement  tends 
to  be  very  eomputationally  expensive  to  later  simulation.  Loeal  grids  refinement  approaeh  ean  also 
be  adapted,  but  extra  interpolation  proeedure  has  to  be  implemented  for  different  levels  of  grids.  By 
generating  enough  grids  points  near  the  irregular  interfaee,  the  order  of  aeeuraey  for  eut-eell  method 
is  maintained.  The  seeond  approaeh  is  to  reduee  the  number  of  grids  in  the  loeal  finite-differenee 
steneil,  where  eorresponding  lower  order  of  non-uniform  finite-differenee  method  ean  be  applied. 
By  using  this  treatment,  the  loeal  order  of  aeeuraey  is  redueed.  If  the  boundary  nearly  degenerates, 
the  aeeuraey  of  eurrent  eut-eell  method  will  not  be  aflfeeted,  provided  that  the  boundary  interfaee 
ean  be  expressed  aeeurately  in  a  numerieal  or  analytieal  form. 
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5,2.2.3.  Boundary  points 


As  shown  in  Fig.  11,  boundary  points  are  the  marker  points  loeated  at  the  interseetion  between  the 
roughness  surfaee  and  grid  lines.  The  governing  equations  are  not  solved  at  these  points.  However, 
flow  variables  at  these  points  are  needed  for  flnite-differenee  algorithms  of  the  irregular  points  near 
the  roughness  surfaee.  The  solid  interfaee  imposes  the  non-slip  and  non-flow-through  eonditions  for 
flow  veloeities  at  the  boundary  points.  Depending  on  the  aetual  thermal  eonditions  of  the  solid 
surfaee,  either  an  isothermal  or  adiabatie  wall  ean  be  speeifled.  In  this  paper,  only  the  isothermal 
wall  is  eonsidered.  In  this  ease,  both  the  veloeities  and  temperature  of  the  boundary  points  are 
speeifled  by  the  wall  boundary  eonditions.  The  loeal  pressure  on  the  solid  surfaee  (boundary  points) 
needs  to  be  eomputed  by  the  loeal  flow  eonditions  near  the  wall. 

There  are  several  approaehes  to  eompute  the  pressure  at  the  boundary  points.  One  approaeh  is  to 
integrate  a  loeal  wall-normal  momentum  equation  to  obtain  the  wall  pressure.  We  ean  also  use  an 
approximation  assumption  of  zero  pressure  gradients  at  the  wall  to  determine  the  pressure  there.  In 
a  previous  paper,  Zhong  used  a  flfth-order  polynomial  extrapolation  to  determine  the  wall  pressure 
[1].  Satisfaetory  results  have  been  obtained  with  this  approaeh.  We  follow  the  extrapolation 
approaeh  of  Zhong  [1]  to  determine  the  pressure  at  the  boundary  points. 

In  eomputing  pressure  at  the  boundary  points,  the  polynomial  extrapolation  is  required  to  have 
eomparable  order  of  aeeuraey  as  that  of  the  interior  sehemes  in  order  to  maintain  the  expeeted 

global  aeeuraey.  To  maintain  a  {p-\-\)-th  order  global  aeeuraey  of  the  upwind  sehemes  developed 

in  the  last  seetions,  it  is  desirable  to  have  at  least  p-th  order  of  aeeuraey  for  the  extrapolation  for 
the  boundary  points. 

Sinee  a  boundary  point  is  formed  by  the  interseetion  of  the  roughness  interfaee  with  one  of  the  grid 
lines,  the  extrapolation  is  eondueted  along  the  direetion  of  the  same  grid  line.  We  use  the  grid  line 

along  the  direetion  as  an  example.  The  methods  ean  be  applied  the  other  direetions  similarly.  In 

two  and  three-dimensional  eases,  there  is  options  of  either  eonstrueting  the  extrapolation  along  the 
direetion  normal  to  the  solid  interfaee,  or  doing  it  along  the  grid  lines.  Though  either  method  ean  be 
used  in  multi-dimensional  problems,  the  latter  approaeh  is  used  in  this  paper  sinee  it  is  simper  in 
implementation  and  more  stable  for  simulations  of  eurrent  tests.  In  this  ease,  the  determination  of 
pressure  at  the  boundary  points  is  a  one-dimensional  extrapolation  along  one  of  the  grid  lines.  The 
one-dimensional  steneil  of  Fig.  13,  whieh  involves  non-uniform  grid  spaeing  9 ,  is  used  to  derive 
the  extrapolation  formulas  for  the  pressure  at  the  boundary  point.  For  example,  for  the  ease  of  four 

point  extrapolation  shown  in  Fig.  13,  pressure  at  the  boundary  point  p{^y)  is  obtained  by  a  third 
degree  polynomial  interpolating  through  the  following  four  interior  pressures:  p{^2)  ^  ^ 

pi^^),  and  pi^s) .  The  grid  spaeing  between  neighboring  grid  points  is  a  eonstant  value  of  Ah  , 
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with  the  exception  that  the  distance  between  the  first  and  second  points  are  6  .  If  there  is  a  dropped 
point  in  the  stencil,  the  dropped  point  is  not  used  in  the  extrapolation  calculations. 


For  two-dimensional  problems,  the  high-order  extrapolation  along  thci^  or  rj  direction  is  employed 

to  extrapolate  pressure  from  the  interior  domain  into  the  boundary.  A  total  of  p  grid  points  in  the 

direction  associated  with  the  boundary  point  are  chosen  to  be  included  in  the  extrapolation  stencil. 
Again,  dropped  points  are  not  included  in  the  extrapolation  stencil.  This  procedure  can  prevent  the 
small  cell  problem  of  producing  numerical  instability  near  the  solid  boundary. 

5,2,2,4,  Dropped  points 

In  this  paper,  in  order  to  avoid  the  small  cell  problem,  a  grid  point  is  designated  as  a  dropped  point 
if  its  distance  from  a  neighboring  boundary  point  along  a  grid  line  is  smaller  than  a  pre-determined 
value  0 .  The  dropped  point,  which  is  associated  with  a  grid  direction,  is  removed  in  the 
finite-difference  grid  stencil  along  the  grid  line  for  the  irregular  points  near  the  boundary.  On  the 
other  hand,  a  grid  point  may  become  a  dropped  point  in  one  direction,  but  remain  a  regular  or 

irregular  point  in  another.  For  example,  point  in  Fig.  14  is  a  dropped  point  in  the  ^  direction, 

but  a  regular  point  in  the  rj  direction.  In  this  case,  the  flow  variables  at  point  FJ  are  not  used  in 

finite  difference  formulas  for  derivatives  in  the  ^  direction.  However,  the  flow  variables  at  the 

same  grid  point  are  needed  for  finite  difference  formulas  for  derivatives  in  the  rj  direction.  The 

flow  variables  at  this  dropped  point  Pj  are  obtained  by  the  interpolation  of  a  stencil  along  the  ^ 
direction. 

For  a  grid  point  which  is  a  regular  or  irregular  point  in  one  direction,  but  a  dropped  point  in 
another,  finite  difference  schemes  along  the  former  direction  may  include  this  point  in  its  stencil. 

As  shown  in  Fig.  14  for  the  case  of  p  =  3,  the  finite-difference  stencil  for  an  irregular  point  Q 
located  at  contains  five  points  in  the  rj  direction,  which  are  in  set 

=  (^,,7,-2)}  •  There  are  two  dropped  points  along  the  ^ 

direction  in  this  set  of  stencil  Qg  :  point  Py  at  and  point  Pj  .  If  the  points 

Pj  and  P2  are  removed  from  the  stencil  used  in  the  calculations  of  flux  terms 
dF  !  dt]  and  dF^  I  dt]  in  g ,  the  stencil  set  Qg  needs  to  be  shifted  two  grids  down  to  include 
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(<^,,77,_3)  and  (<^,,?7,_4)  to  maintain  the  accuracy.  In  this  case,  the  resulting  stencil  for  Q  may 

contain  a  significantly  large  interval  6  compared  with  the  normal  grid  spacing  A/z ,  which  may 
lead  to  a  deterioration  of  accuracy  of  the  method.  Therefore,  we  maintain  the  original  grid  stencil 

Qg ,  which  includes  points  Pj  and  ,  along  the  77  direction  for  point  Q  .  We  calculate  the  flow 
variables  of  these  two  dropped  points  by  interpolation  along  the  ^  direction. 
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Fig.  14.  An  example  of  a  grid  point,  such  as  point  TJ  and  point  ,  which  is  a  regular  or  irregular  point  in  one 

direction,  but  a  dropped  point  in  another  for  the  case  of  p  =  3'.  V  represents  regular  points,  0  dropped 

points,  O  boundary  points,  and  •  irregular  points. 

For  the  case  of  p  =  3  for  point  Pj ,  a  third-order  polynomial  interpolation  along  the  ^  direction 
is  employed  to  compute  the  flow  variables  at  this  point.  The  interpolation  is  carried  out  along  the 
direction,  which  is  along  the  direction  where  the  point  is  dropped.  For  p  =  3  as  shown  in  Fig. 

14,  the  interpolation  stencil  for  point  P,  is  set  Qp  =  B}  ,  where  B 

represents  the  boundary  point.  For  higher  order  method,  the  order  of  interpolation  needs  to  be 
increased  accordingly.  For  p-4,  a  fourth  order  interpolation  should  be  used.  For  a  general  case 

of  p-th  order  methods  at  the  boundary,  a  total  of  p-\  adjacent  grid  points  and  exactly  one 

boundary  point  along  the  direction  are  chosen  as  the  interpolation  stencil.  The  interpolant  can 
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be  written  as 


z 


«=1 


f  /r  4zil 


U„ 


(19) 


Where  is  the  eoordinate  of  dropped  point,  {U.,i  =  \..p]  is  eonservative  flow  variables  at 
the  grid  points  of  the  interpolation  steneil. 


Similar  interpolation  proeedures  ean  also  be  earried  out  if  a  point  is  designated  as  a  dropped  point 
in  the  J]  direetion,  but  is  an  irregular  or  regular  point  in  the  ^  direetion.  If  a  grid  point  is 

designated  as  a  dropped  point  in  both  the  and  r/  direetions,  there  is  no  need  to  do 
interpolations  beeause  this  point  is  removed  from  the  ealeulations  of  both  direetions. 


5,3  Test  and  applications 

The  new  high-order  eut-eell  method  has  been  tested  in  the  eomputations  of  several  one  and 
two-dimensional  hyperbolie  and  elliptie  equation  in  irregular  domains.  The  results  show  that  up  to 
fourth-order  aeeuraey  in  both  and  norm  ean  be  obtained  for  the  eurrent  eut-eell  method  for 
both  problems.  We  subsequently  have  applied  a  third-order  eut-eell  method  to  the  two-dimensional 
eompressible  Navier-Stokes  equations  for  simulating  roughness  indueed  reeeptivity  for  hypersonie 
flow  over  a  flat  plate  with  a  blowing  and  suetion  slot  plaeed  near  the  leading  edge.  The  surfaee 
roughness  height  is  approximately  half  of  the  loeal  boundary  thiekness.  By  using  the  uniformly 
third-order  eut-eell  method,  small  disturbanees  generated  by  the  blowing  and  suetion  slot  are  well 
resolved  around  the  roughness  element.  The  results  obtained  by  the  new  third-order  eut-eell  method 
are  eonsistent  with  those  of  the  linear  stability  analysis  results.  In  addition,  steady  solutions  of  the 
flow  aflfeeted  by  the  roughness  are  well  eaptured  by  the  eut-eell  method.  The  simulation  results 

suggest  that  the  existenee  of  small  surfaee  roughness  (  )  affeets  the  propagation  of 

instability  waves  in  the  hypersonie  boundary  layer.  The  details  regarding  the  tests  ean  be  fond  in  [1]. 


6  STABILIZATION  OF  A  MACH  5.92  FLAT-PLATE  BOUNDARY  LAYER 
WITH  FINITE  HEIGHT  SURFACE  ROUGHNESS 


One  of  the  major  aeeomplishments  of  the  eurrent  researeh  projeet  is  on  revealing  the  reeeptivity 
meehanism  of  a  hypersonie  boundary  layer  with  finite  surfaee  roughness  to  wall  blowing-suetion. 
These  results  are  diseussed  in  more  detail  in  this  seetion.  Details  of  these  results  have  been 
published  in  Publications:  [3],  [9],  [13],  [15],  [20], 
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6,1  Flow  conditions  and  roughness  model 


With  the  development  of  advaneed  eomputers  and  numerieal  teehniques,  numerieal  simulation  of 
the  reeeptivity  proeess  by  direetly  solving  Navier-Stokes  equations  has  beeome  feasible.  By  solving 
the  eompressible  linearized  Navier-Stokes  equations,  Malik  et  al.  [65]  investigated  the  responses  of 
a  Maeh  8  flow  over  a  sharp  wedge  of  a  half-angle  of  5.3°  to  three  types  of  external  foreing:  a  planar 
free-stream  aeoustie  wave,  a  narrow  aeoustie  beam  enforeed  on  the  bow  shoek  near  the  leading 
edge,  and  a  blowing-suetion  slot  on  the  wedge  surfaee.  They  eoneluded  that  these  three  types  of 
foreing  eventually  resulted  in  the  same  type  of  instability  waves  in  the  boundary  layer.  However  the 
reeeptivity  meehanism  was  not  studied  in  detail.  Ma  and  Zhong  [66]  studied  the  reeeptivity 
meehanisms  of  the  same  hypersonie  boundary  layer  to  various  freestream  disturbanees,  i.e.,  fast  and 
slow  aeoustie  waves,  vortieity  waves,  and  entropy  waves,  by  solving  the  two-dimensional 
eompressible  Navier-Stokes  equations.  They  found  that  the  stable  modes  in  the  boundary  layer 
played  a  very  important  role  in  the  reeeptivity  proeess. 

In  this  report,  the  new  eut-eell  method  introdueed  in  Seetion  5  is  used  to  eompute  hypersonie 
viseous  flow  over  a  flat  plate  with  an  isolated  surfaee  roughness  element  (Fig.  1).  Both  steady  and 
unsteady  flows  are  eonsidered.  The  freestream  flow  eonditions  are  the  same  as  those  used  in 
Maslov’s  experiment  [13]  as  follows, 

jp,=0.72,i?^=p^*u„///^=1.32xl0Vm 
=  5.92,r„  =  48.69i:,P^  =  742.76Pa 

where  ,  T^,  P^,  P^,  are  Maeh  number,  temperature,  pressure,  Prantle  number  and  unit 

Reynolds  number,  respeetively.  The  flat  plate  is  assumed  to  be  isothermal  with  a  eonstant 
temperature  of  =  350. 07f  .  The  total  length  of  the  flat  plate  is  about  1 .69m  . 


In  this  seetion,  unless  stated  otherwise,  most  flow  variables  are  presented  as  dimensionless  ones. 
The  main  exeeption  is  that  the  x  and  y  eoordinates  are  presented  in  dimensional  form  so  that  the 
simulational  results  ean  be  easily  related  to  the  experimental  setup.  For  other  variables,  we 
nondimensionalize  the  flow  veloeities  by  the  freestream  veloeity  ,  density  by  ,  pressure  by 

Poo  (^00  >  and  temperature  by  . 

An  isolated  roughness  element  of  smooth  shape  is  plaeed  on  the  surfaee  of  the  flat  plate  at 
x  =  0.185w2  downstream.  Motivatied  by  Whiteheard’s  experiments  [67],  the  shape  of  the  surfaee 
roughness  is  ehosen  to  be  a  two-dimensional  bump,  governed  by  the  following  elliptie  equation. 
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The  computations  of  this  study  are  performed  under  the  following  parameters: 
a  =  2,  Z)  =  1,  =  0.185m  and  h  =  5l2  where  =0. 00081m  which  corresponds  to  the  boundary 

layer  thickness  at  x^  as  shown  in  Fig.  15. 


Fig.  15.  1)  Computational  grid  for  hypersonic  flow  over  a  flat  plate  with  an  isolated  surface  roughness.  The 
grid  size  is  241x121  (showing  one  out  of  four  grid  points  in  both  directions,  2).  Comparison  of  the  current 

shock-fitting  results  with  theoretical  solutions  {T^  and  u^)  for  temperature  and  wall  normal  velocity 

distributions  long  a  wall-normal  grid  line  located  at  x  =  0.1676m  on  the  plate  surface. 


In  this  report,  a  third-order  cut-cell  method  described  in  Section  5  is  used  to  compute  the 
two-dimensional  viscous  hypersonic  flow  over  the  flat  plate  with  the  roughness  element.  As 
described  in  Section  5,  a  coordinate  transformation  is  employed  to  transform  the  physical  domain 
shown  in  Fig.  15  into  a  rectangle  computational  domain  with  a  set  of  Cartesian  grid.  The  optimal 
transformation  formula  is  determined  by  the  specific  physical  problem  considered.  For  viscous  flow 
over  a  flat  plate,  it  is  necessary  to  cluster  more  grid  points  at  the  bottom  wall  surface  in  order  to 
resolve  the  viscous  boundary  layer.  In  this  paper,  a  two-step  mapping  procedure  is  used  to  obtain 
better  resolution  inside  the  viscous  boundary  layer. 

In  the  first  step  of  the  transformation,  the  entire  physical  domain  shown  in  Fig.  15  is  transformed 
into  a  square  domain  defined  on  [0,l]x[0,l]  as  an  intermediate  coordinate  space.  The 
transformation  relation  is  defined  as: 
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where  (X,Y)  is  defined  under  the  intermediate  eoordinate  system,  L  is  the  streamwise  length  of 

the  flat  plate  in  physieal  domain,  is  the  streamwise  eoordinate  of  starting  point  of  eomputation 

in  eaeh  zone,  and  the  distanee  along  r/  direetion  between  the  solid  wall  and  bow  shoek  is  referred 
to  H. 

In  the  seeond  step  of  the  transformation,  the  intermediate  plane  in  the  (X,  Y)  spaee  is  mapped  into 

the  final  eomputational  domain  in  order  to  eluster  more  grid  points  into  the  viseous 

boundary  layer  near  the  flat  plate  surfaee.  In  the  present  study,  an  exponential  streteh  funetion  is 
used  to  eluster  grids  in  the  t]  direetion  as  follows. 


r  = 
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(23) 


where  P  is  the  stretehing  parameter.  The  value  of  y^  =  1.01is  used  in  this  paper.  With  this 

P  value,  about  50%  of  the  total  grid  points  in  the  rj  direetion  are  elustered  inside  the  boundary 

layer.  The  eombination  of  the  transformations  of  these  two  steps  leads  to  overall  transformation 
given  by  Eq.  (2). 

In  the  eomputational  domain  (p,  t]),  the  baseline  grid  is  a  set  of  Cartesian  grid  similar  to  the 

sehematie  shown  in  Fig.  1 1  (2).  The  interaetion  of  the  roughness  surfaee  with  the  grid  lines  ereates 
eut  eells,  whieh  are  treated  by  the  eurrent  eut-eell  method.  The  surfaee  equation  (21)  of  the 
roughness  surfaee  is  transformed  into  the  eomputational  domain  in  order  to  apply  the  eut-eell 
method.  By  substituting  Eqs.  (22)  and  (23)  into  the  roughness  surfaee  Eq.  (21),  the  analytieal 
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equation  for  roughness  surfaee  in  the  eomputational  domain  ean  be  written  as 


f(P  I  m^L+Xs..f 

^2  +  jj2 


where  f  = 
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From  Eq.  (24),  we  ean  ealeulate  the  eoordinates  of  the  boundary  points  by  eomputing  the 
eoordinates  of  the  interseetion  between  the  roughness  surfaee  and  the  grid  lines.  For  a  given  value 

of  along  the  roughness  surfaee,  the  Newton  iteration  method  is  employed  to  solve  Eq.  (24)  for 

the  eorresponding  value  of  rj  on  the  surfaee  for  a  boundary  point,  and  viee  versa.  Having 

eomputed  the  boundary  points,  we  ean  identify  all  irregular  and  dropped  points  near  the  roughness 
surfaee.  The  remaining  grid  points  are  regular  points. 

After  eomputing  all  eoordinates  of  the  irregular,  regular,  boundary  and  dropped  points  as  shown  in 
Fig.  11,  the  eurrent  high-order  eut-eell  method  is  applied  to  diseretize  Eq.  (3).  The  bow  shook  is 
treated  by  a  high-order  shook- fitting  method.  The  overall  aoouraoy  of  simulation  is  determined  by  the 
order  of  soheme  used  for  eomputing  the  invisoid  and  visoous  flux  terms. 

6,2.  Steady  flow  solution  without  surface  roughness 

The  steady  mean  flow  solutions  are  oaloulated  by  using  the  fifth-order  shook-fitting  method.  But  in 
the  leading  edge  region  of  the  fiat  plate,  there  is  a  singularity  at  the  tip  of  the  plate  and  the 
high-order  shook  fitting  method  oannot  be  used  there.  Thus  a  seoond-order  TVD  shook-oapturing 
method  is  employed  to  ealeulate  a  small  looal  flow  field  around  the  tip  of  the  flat  plate.  The 
eomputational  domain  for  the  TVD  oaloulations  starts  at  x  =  -0.006m  and  ends  at  a  very  short 
distanoe  downstream  of  the  leading  edge  at  x  =  0.003m  .  A  total  241x121  grid  points  are  used.  A 
TVD  soheme,  whioh  follows  that  used  by  Lee  et  al.  [68],  is  applied  to  Eq.  (3).  The  semi-disorete 
system  of  ordinary  differential  equations  are  then  solved  by  using  a  fourth-order  Runge-Kutta 
method. 

Having  obtained  the  steady  state  solutions  at  the  leading  edge,  we  then  use  the  solution  of  the  TVD 
soheme  as  the  inlet  oondition  to  start  the  subsequent  shook-fitting  oaloulations.  The  eomputational 
domain  for  the  high-order  shook-fitting  methods  starts  at  x  =  0.003m  and  ends  at  x  =  1.68784m  .  In 
aotual  simulations,  the  eomputational  domain  is  divided  into  30  zones,  with  total  of  5936  grid  points 
in  the  streamwise  direotion  and  121  points  in  the  wall-normal  direotion.  As  mentioned  above,  the 
seoond  zone  uses  the  results  of  the  first  zone  of  the  seoond-order  TVD  solution  as  the  inlet  oondition. 
A  later  zone  uses  the  interpolation  of  its  former  zone’s  data  as  the  inlet  oondition. 
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The  solutions  for  the  steady  viseous  flow  over  the  flat  plate  with  no  surfaee  roughness  are  first 
obtained.  Fig.  15  shows  the  streamwise  veloeity  and  temperature  profiles  along  the  wall-normal 
direetion  at  x  =  0.1676m  near  the  wall.  The  eurrent  numerieal  solutions  obtained  by  the  fifth-order 
shoek-fitting  seheme  are  eompared  with  the  self-similar  boundary  layer  solution.  In  order  to 

eompare  with  the  self-similar  solution  of  the  boundary  layer,  the  y  eoordinate  is 

nondimensionalized  by  .  Fig.  15  illustrates  that  the  results  of  the  eurrent  numerieal 

simulation  agree  very  well  with  the  theoretieal  solutions  near  the  wall.  The  seeond-order  TVD 
seheme  is  aeeurate  enough  to  be  used  as  the  inlet  eondition  of  the  fifth-order  shoek-fitting  method 
downstream  of  the  leading  edge. 

6,3.  Steady  flow  solution  with  an  isolated  surface  roughness 

The  roughness  model  proposed  in  See  6.1  is  simulated  by  using  the  high-order  eut-eell  method  with 
same  flow  eondition  in  See  6.2.  In  this  report,  the  results  of  a  third-order  {p  =  !>)  eut-eell  method  is 
presented.  In  simulation,  the  numerieal  dissipation  eoeffieient  a  is  ehosen  within  range 
[0.25  -  6.0] .  All  the  simulations  are  very  stable.  Thus  the  numerieal  dissipation  of  the  interior  seheme 
may  only  have  minor  influenee  on  the  overall  stability. 

Fig.  16  shows  steady  flow  solutions  eomputed  by  the  third-order  eut-eell  method,  in  both  the 
streamline  pattern  and  the  eontours  of  wall  normal  veloeity  and  pressure  eomponents.  In  order  to 
assess  the  numerieal  aeeuraey  of  the  solutions,  steady  solutions  are  obtained  by  using  the  following 
two  sets  of  grids:  241x121  and  441x121.  Due  to  the  restrietion  of  eomputational  time,  the  grid 
refinement  study  is  only  performed  with  respeet  to  streamwise  resolution  in  one  zone.  The 
wall-normal  veloeities  eomputed  by  fine  and  eoarse  grids  are  plotted  with  the  same  levels  of 
eontour  line  in  Fig.  16.  The  good  agreement  between  these  two  sets  of  results  suggests  that  the  grid 
resolution  used  in  the  eurrent  simulations  is  high  enough  for  the  steady  flow  simulation.  Fig.  16  also 
shows  the  steady-flow  pressure  eontours  for  flow  over  a  flat  plate  with  an  isolated  roughness 
eomputed  by  the  third-order  eut-eell  method.  Beeause  the  flow  is  supersonie  behind  the  bow  shook, 
a  family  of  Maoh  waves  is  generated  by  flow  over  the  roughness.  The  oompression  waves  are 
followed  by  expansion  waves  when  the  flow  expands  around  the  roughness  surfaee.  These  Maoh 
waves  are  approximately  parallel  to  the  bow  shook  interfaoe  in  the  later  zones. 
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Fig.  16.  Hypersonic  flow  soiutions  in  the  region  near  the  roughness  eiement  computed  by  the  third-order 
cut-ceii  method:  1)  streamiine  pattern,  2)  contour  iines  of  waii  normai  veiocity  components  computed  by  using 
two  sets  of  grids  (coarse  grid:  241x121 ,  fine  grids:  441x121 ),  3)  waii  normai  veiocity  contour  4), Pressure 
contours  for  two-dimensionai  hypersonic  flow  over  a  fiat  piate  with  an  isoiated  roughness  computed  by  the 
third-order  cut-ceii  method.  The  totai  grid  points  in  this  range  of  the  flow  fieid  is  1920x121 . 


The  steady  flow  solutions  demonstrate  that,  affeeted  by  the  roughness,  hypersonie  boundary  layer 
over  the  flat  plate  is  modified  significantly  in  the  regions  both  downstream  and  upstream  of  the 
roughness  element.  The  parallel  flow  assumption  is  no  longer  valid  for  flow  near  the  roughness. 
Consequently,  normal-mode  linear  stability  analysis  may  not  be  accurate  in  this  flow  region.  Thus  a 
different  boundary  layer  instability  mechanism  may  be  introduced.  One  the  other  hand,  the  effects 
of  the  roughness  is  most  significant  in  the  regions  surrounding  the  roughness.  The  roughness 
effects  on  the  steady  flow  decay  if  we  move  further  downstream.  Thus  the  flow  is  approximately 
parallel  again  in  the  far  downstream  of  the  surface  roughness.  Linear  Stability  Theory  (LST)  can  be 
applied  to  analyze  the  flow  field  which  is  approximately  parallel.  Good  agreement  between 
numerical  and  theoretical  solutions  can  be  achieved.  More  details  can  be  found  in  [69]. 


6,4,  Receptivity  of  hypersonic  flow  with  surface  roughness  to  wall  blowing  and  suction 

In  this  section,  the  new  third-order  cut-cell  method  is  applied  to  the  computation  of  transient 
responses  of  the  same  Mach  5.92  boundary  layer  to  forcing  waves  introduced  by  a  blowing  and 
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suction  slot  located  on  the  plate  surface.  It  is  termed  the  receptivity  problem  of  the  boundary  layer 
to  forcing  waves  [24],  The  reeeptivity  problem  is  eritical  to  the  understanding  of  physical 
mechanisms  of  hypersonic  boundary  layer  transition.  The  reeeptivity  of  the  same  hypersonic 
boundary  layer  to  wall  blowing  and  suction  has  been  studied  by  Wang  and  Zhong  [70]  and  [71]  for 
smooth  surfaee  without  roughness.  In  this  paper,  we  study  the  additional  effects  of  the  isolated 
surface  roughness  on  the  reeeptivity  process  by  using  the  eurrent  third-order  eut-cell  method.  The 
current  results  are  compared  with  those  of  [70]  and  [71]  to  study  the  effeets  of  surface  roughness  on 
the  reeeptivity  process. 

The  reeeptivity  study  is  mainly  eoneerned  with  the  excitation  of  instability  waves,  the 
characteristics  of  which  can  be  analyzed  by  the  linear  stability  theory  [72].  The  instability  theory 
analyzes  the  propagation  of  individual  sinusoidal  waves  in  the  streamwise  direction  inside  the 
boundary  layer.  These  waves  are  referred  as  Tollmien-Schliehting  (T-S)  waves  for  low  speed  flow. 
The  instability  waves  are  vortieity  waves,  whose  amplitudes  vary  though  the  boundary  layer  and  die 
off  exponentially  outside  the  boundary  layer.  For  the  case  of  small  perturbations  in  the  flow  field, 
the  perturbations  of  flow  variables  can  be  written  in  the  form  of  a  normal  mode,  i.e.. 
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where  6  is  the  perturbation  of  temperature.  The  frequeney  of  the  waves  is  co  and  the  streamwise 
and  spanwise  wave  numbers  of  the  instability  wave  are  represented  by  a  ,  and  /?  respectively. 


Substituting  Eq.  (25)  into  a  linearized  version  of  the  full  Navier-Stokes,  we  obtain  a  system  of 
linearized  equations  of  the  stability  theory.  These  equations  are  required  to  satisfy  a  number  of 
boundary  eonditions.  By  imposing  the  homogenous  physieal  eonditions,  the  number  of  solutions  of 

an  eigen  problem  to  the  linearized  equations  with  specifie  value  of  a,  (3  and  W  is  eonstrained. 


The  relation  for  the  instability  wave  parameter  a,  (3  and  W  are  referred  as  dispersion  relations  in 
the  following  form, 

W  =  W{a,l3)  (26) 

Extensive  numerieal  and  theoretieal  researeh  has  been  eondueted  to  solve  the  linearized 
Navier-Stokes  equations  and  many  eharaoteristies  regarding  the  instability  waves  in  hypersonie 
boundary  layers  have  been  diseovered.  [13,  72-75].  Mack  [72]  identified  the  unstable  modes  by 
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using  the  linear  stability  theory.  He  showed  that  inside  a  supersonie  boundary  layer,  there  are 
multiple  higher  instability  modes  in  addition  to  the  first  mode,  whieh  is  the  eompressible 
eounterpart  of  T-S  waves  in  the  ineompressible  boundary  layers.  These  instability  modes  in  the 
supersonie  boundary  layer  are  termed  as  first  mode,  seeond  mode,  third  mode,  ete.  The  seeond 
mode  is  also  ealled  the  Maek  mode.  For  supersonie  boundary  layer  with  Maeh  number  larger  than 
four,  the  seeond  Maek  mode  is  the  most  unstable  mode,  and  it  plays  an  important  role  in 
hypersonie  boundary  layer  transition. 

To  exeite  the  propagation  of  small  disturbanees  inside  the  boundary  layer,  a  blowing  and  suetion 
slot  is  imposed  as  periodie-in-time  boundary  eonditions  for  the  perturbations  of  the  mass  flux  on  the 
wall.  The  blow-suetion  slot  is  loeated  at  x  =  0.030m  and  spreads  over  several  grids  spaees.  The 
perturbations  in  the  blowing  and  suetion  slot  are  governed  by  the  following  funetion: 

Pv  =  ^og(/)Isin(w„0  (27) 


where  is  an  amplitude  parameter,  -  nf^  is  eireular  frequeney  of  this  multi-frequeney 
perturbation,  and  g(/)  is  a  non-dimensional  x-direetion  profile  funetion  defined  as 


gil) 


J2O.24/'  -35.4375^+15.1875/'  (/  <  1) 
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(28) 


The  variable  I  is  the  non-dimensional  parameter  assoeiated  with  eurrent  eoordinate  of  this  blowing 
and  suetion. 


(29) 

he  basie  frequeney  fy  is  given  by  f  =  50kHz ,  and  the  other  14  different  frequeneies  are  imposed 

as  /  =  100,150,200...  kf/z  for  n  =  2,3,4....  This  steady  and  unsteady  flow  eonditions  used  in  the 

eurrent  simulation  are  the  same  as  those  used  in  the  simulation  eondueted  by  Wang  et  al.  [70]  for 
the  eases  of  smooth  wall  without  roughness.  Balakumar  [76]  also  investigated  the  reeeptivity  of  a 
2-D  roughness  to  aeoustie  waves  and  found  the  isolated  roughness  does  not  eontribute  mueh  in 
generating  unstable  disturbanees. 
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Fig.  17.  Comparison  of  contours  of  pressure  perturbations  for  unsteady  Mach  5.92  flow  over  a  fiat  piate  with 
and  without  surface  roughness.  The  forcing  frequency  is  f  =  l50kHz  .  The  current  resuits  for  flow  without  1), 
or  with  2)  roughness  are  computed  by  the  third-order  cut-ceii  method. 


The  simulation  is  carried  out  with  combined  forcing  waves  of  fifteen  frequencies.  A  Fast  Fourier 
Transformation  (FFT)  technique  is  used  to  separate  the  results  with  different  frequencies.  Fig.  17 
compares  the  contours  of  pressure  perturbations  for  unsteady  Mach  5.92  flow  over  a  fiat  plate  with 

and  without  surface  roughness.  The  forcing  frequency  is  /  =  1 50kHz  .  The  current  results  for  flow 

with  roughness  are  computed  by  the  third-order  cut-cell  method.  The  corresponding  results  for  the 
smooth  surface  without  roughness  was  obtained  by  Wang  and  Zhong  [70]  and  [71].  Fig.  17  (1) 
shows  that,  for  the  case  of  no  surface  roughness,  instability  waves,  which  are  induced  by  the  forcing 
waves  introduced  at  the  blowing  and  suction  slot,  propagate  dowstream  in  the  streamwise  direction 
inside  the  boundary  layer.  Fig.  17  (2)  shows  how  the  roughness  element  interacts  with  the  wave 
propagation  inside  the  boundary  layer. 


Fig.  18  compares  the  contours  of  pressure  perturbations  for  cases  with  and  without  surface 
roughness  for  a  longer  flow  region  with  combined  15  forcing  frequencies.  The  roughness  is  located 

at  =  0.185m .  The  pressure  disturbances  at  flat  plate  surface  are  reduced  significantly  after  the 

instability  waves  pass  the  surface  of  roughness  element.  Thus  the  local  surface  roughness  could 
potentially  stabilize  the  instability  waves  in  the  boundary  layer. 
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Fig.  18.  Comparison  of  contours  of  pressure  perturbations  for  unsteady  Mach  5.92  flow  over  a  fiat  piate  with 
and  without  surface  roughness  for  a  ionger  flow  region.  The  forcing  frequency  is  combined  15  frequencies. 

The  roughness  is  iocated  at  =  0.185m  . 

A  muti-domain  spectral  collocation  (MDSC)  method  of  Malik  [77]  is  used  to  conduct  stability 
analysis  of  the  steady  solution  for  the  hypersonic  flow  over  the  roughness  element.  More  details 
about  the  linear  stability  analysis  and  methods  can  be  found  in  [77]. 

The  dimensionless  frequency  used  for  linear  stability  analysis  is  defined  as 

F  =  (30) 

«oo 


where  F  is  dimensionless  frequency,  v  is  kinematic  viscous  coefficient.  In  present 
simulation,^  =  6.05  X 10  m  Is^  F  =  5. 30x10  for  the  case  of  and  ^co  -  827.29m/ 5 

In  LST  analyses  of  boundary  layer  flows,  the  Reynolds  number  based  on  the  focal  length  scale  of 
boundary  layer  thickness,  5  ,  is  generally  used.  They  are  expressed  as 


_  PJ^J  I 


(31) 


Hence  the  relation  between  the  R  and  the  unite  Reynolds  number  is. 


R  =  y[R^ 


(32) 


With  the  definitions  of  Reynolds  number  R  and  the  dimensionless  frequency  F ,  the  dimensionless 
circular  frequency  m  can  also  be  expressed  as 


co  =  RF 


(33) 
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Fig.  19.  Wave  mode  profiles  obtained  by  linear  stability  analysis  at  various  locations:  1)  mode  F  profile  for 
velocity,  2)  mode  F  profile  for  pressure,  3)  mode  S  profile  for  velocity,  4)  mode  S  profile  for  pressure.  The 
amplitudes  of  disturbances  are  nondimensionlized  by  their  corresponding  values  on  the  surface  of  flat  plate. 

There  are  two  major  wave  modes  in  hypersonie  boundary  layers:  mode  F  and  mode  S.  Fig.  19 
shows  wave  mode  profiles  obtained  by  linear  stability  analysis  at  various  loeations  for  both  mode  F 
and  mode  S.  The  disturbanee  amplitudes  are  nondimensionlized  by  the  pressure  disturbanee  on  the 

surfaee  of  flat  plate,  e.g.  u  '(y)  =  u{y)  /  ^(0) .  In  the  upstream  of  roughness  at  x  =  0. 1620/n  where  the 

roughness  effeets  on  the  mean  flow  are  not  that  signifleant,  mode  S  and  mode  F  ean  be  indentifled 
elearly.  As  we  move  gradually  downstream,  the  profiles  of  the  wave  modes  vary  signifieantly. 
After  X  =  0.1980m ,  mode  S  and  mode  F  ean  be  identified  again  by  examining  their  perturbation 
profiles.  The  mean  flow  solution  obtained  by  the  eut-eell  method  ean  produee  aeeuraey  results  for 
the  stability  analysis. 

6,5  Stabilization  effect  of  location  of  the  finite  roughness 

In  order  to  investigate  the  effeet  of  loeation  of  the  finite  surfaee  roughness  on  reeeptivity  and 
stabilization,  a  series  of  numerieal  simulations  have  been  earried  out  for  different  roughness 
loeations.  Speeifieally,  four  oases  of  different  roughness  loeations  are  oonsidered.  Fig.  20  shows  a 
sehematie  of  different  roughness  loeations  for  the  four  eases.  In  eaeh  ease,  wall  blowing-suetion 
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with  fifteen  frequeneies  is  introdueed  on  the  flat  plate  surfaee.  The  subsequent  responses  of  the 
boundary  layer  are  simulated  by  the  high-order  eut-eell  finite  differenee  method. 
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Fig.  20  A  schematic  of  iocations  of  the  biowing-suction  actuator  for  the  four  cases:  1 .  fiat  piate;  2.  bow  shock; 

3.  surface  roughness;  4.  boundary  iayer;  5.  biowing  suction. 


Specifically,  roughness  locations  with  respect  to  the  synchronization  point  at  x  =  0.33  Im  are  as 
follows, 

•  Case  1:  Upstream  of  the  synchronization  point  x  =  0.1 85m  ;  The  roughness  height 

/z  =  i^  =  0.00081m. 

2 

•  Case  2:  Upstream  of  the  synchronization  point  x  =  0.260m  ;  The  roughness  height 

/z  =  1^  =  0.001 11m. 

2 


•  Case  3: 


Synchronization  point  x  =  0.33  Im ;  The  roughness  height 


/2  =  i^  =  0.00141m 
2 


•  Case  4:  Downstream  of  the  synchronization  point  x  =  0.410m;  The  roughness  height 

/2  =  i^  =  0.00171m. 

2 

Fig.  21  shows  the  pressure  amplitude  on  the  flat  plate  surface  with  roughness  located  in  different 
locations.  It  shows  that  for  all  four  cases,  the  development  of  two-dimensional  wall  perturbations 
eventually  result  in  the  same  type  of  instability  wave  (mode  S)  in  the  boundary  layer.  When  the 
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surface  roughness  located  in  the  far  upstream  of  the  synchronization  point  (Case  1),  the  instability 
wave  develops  in  the  same  way  as  the  case  without  roughness.  But  as  the  roughness  moves 
gradually  into  the  downstream  (Case  2),  the  amplitude  of  pressure  in  the  flat  plate  surfaee  is 
amplified.  Mode  S  is  destabilized  when  the  roughness  element  is  plaeed  suffieient  elose  to  the 
synehronization  point.  In  Case  3  and  Case  4,  when  the  roughness  loeated  in  downstream  of 
synehronization  point,  the  finite  surfaee  roughness  stabilizes  the  unstable  wave  inside  the  boundary 
layer.  The  amplitude  of  pressure  perturbation  is  redueed  in  further  downstream. 


Fig.  21.  Numerical  simulation  of  receptivity  of  the  Mach  5.92  flow  over  a  flat  plate  to  two-dimensional  wall 
perturbations  with  surface  roughness:  Amplitude  comparison  of  pressure  for  four  cases. 


Fig.  22.  Numerical  simulation  of  receptivity  of  the  Mach  5.92  flow  over  a  flat  plate  to  two-dimensional  wall 
perturbations  with  surface  roughness:  comparison  of  numerical  and  LST  growth  rates. 
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Fig.  22  shows  the  comparison  of  numerical  and  LST  growth  rate  a,  as  variation  of  dimensionless 

circular  frequency  o.  The  DNS  result  and  LST  does  not  agree  initially  due  to  the  existence  of 

various  eigen  modes  perturbations.  After® >0.12  (x> 0.33m),  mode  S  becomes  dominant  among 

the  eigen  modes  spectral.  Thus  after  the  Mode  S  fully  develops  in  the  boundary  layer,  the 

theoretical  growth  rate  a,  of  Mode  S  plotted  in  red  lines  agrees  well  with  the  DNS  result.  In  the 

region  surrounding  the  surface  roughness  element,  the  pressure  perturbation  on  the  plate  surface 

oscillates  strongly  for  all  four  cases,  which  drive  the  growth  rate  a,  to  be  infinite.  But  in  the 

downstream  of  surface  roughness,  the  mean  flow  is  restored  and  the  Mode  S  growth  rate  can  be 
calculated.  For  Case  1,  the  growth  rate  does  not  change  much  with  the  roughness  effect.  This  shows 
the  roughness  only  affect  local  mean  flow  and  local  mode  S  instability.  But  for  Case  2,  mode  S  is 
destabilized  before  the  synchronization  point  (®»0.12),  and  develops  in  a  similar  way  to  case 
without  roughness  after  that.  Thus  the  modified  local  mean  flow  amplifies  the  amplitude  of 
perturbation.  One  explanation  to  this  is  the  non-parallel  feature  of  local  mean  flow  with  positive 
wall-normal  velocity.  For  Case  3  and  4,  being  very  different  from  Case  1  and  2,  mode  S  is  stabilized 
after  the  synchronization  point.  This  once  again  can  be  explained  as  the  non-parallel  feature  of  the 
local  mean  flow  with  roughness  in  the  upstream  and  negative  wall-normal  velocity.  More  accurate 
theoretical  analysis  tools  beyond  LST  should  be  adopted  to  analyze  this  non-parallel  region. 

To  summarize,  the  current  simulation  results  have  indicated  that  the  roughness  location  plays  an 
important  role  in  the  developments  of  mode  S  by  the  blowing-suction  actuator.  Mode  S  is 
destabilized  only  when  the  roughness  element  is  located  close  upstream  of  the  synchronization  point. 
On  the  other  hand,  when  the  roughness  element  in  is  downstream  of  the  synchronization  point, 
mode  S  is  stabilized.  This  happens  even  when  the  roughness  is  still  within  the  unstable  region  of 
mode  S.  Therefore,  the  synchronization  point  and  finite  roughness  location  is  critical  to  the 
receptivity  process.  The  relationship  between  the  location  of  the  roughness  and  the  synchronization 
point  suggests  that,  in  order  to  control  or  delay  the  laminar-turbulent  transition  more  efficiently,  the 
roughness  element  should  be  placed  upstream  of  the  synchronization  point  between  mode  S  and 
mode  F. 
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